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■ Abstract 

Using methods from soft-collinear and heavy-quark effective theory, a systematic factor- 
ization analysis is performed for the B — > X s j photon spectrum in the endpoint region 
mfe — 2E 7 = O(Aqcd)- It is proposed that, to all orders in 1/my,, the spectrum obeys 
a novel factorization formula, which besides terms with the structure H J ® S familiar 
from inclusive B — > X u l v decay distributions contains "resolved photon" contributions 
of the form H J <2> S (8> J and Here S and J are new soft and jet 

functions, whose form is derived. These contributions arise whenever the photon couples 
to light partons instead of coupling directly to the effective weak interaction. The new 
contributions appear first at order l/m& and are related to operators other than Qj^ 
in the effective weak Hamiltonian. They give rise to non- vanishing 1/mb corrections 
to the total decay rate, which cannot be described using a local operator product ex- 
pansion. A systematic analysis of these effects is performed at tree level in hard and 
hard-collinear interactions. The resulting uncertainty on the decay rate defined with a 
cut Ery > 1.6 GeV is estimated to be approximately ±5%. It could be reduced by an 
improved measurement of the isospin asymmetry Ao- to the level of ±4%. We see no 
possibility to reduce this uncertainty further using reliable theoretical methods. 



1 Introduction and outline 



The radiative decay B — > X s 7 plays an important role in testing the Standard Model and 
constraining its possible extensions at or beyond the TeV scale. Comparing the predictions 
for the branching ratio of this decay obtained in extensions of the Standard Model with 
experiment provides powerful constraints on the parameter space of many new-physics models 
(see e.g. [1-3] for analyses in the context of the MSSM, and [4] for an overview of several other 
models). The calculation of the B — > X s ^ branching ratio in the Standard Model has been 
pushed to the next-to-next-to-leading order in renormalization-group improved perturbation 
theory [5], leading to the prediction Br(5 X s ~f) = (3.15±0.23) • 1(T 4 for a cut E 7 > 1.6 GeV 
on the photon energy measured in the 5-meson rest frame. A dedicated analysis of cut-related 
effects and uncertainties gives the slightly lower value Br(B — > X^) = (2.98 ± 0.26) ■ 10 -4 
[6]. These theoretical estimates are in good agreement with the current experimental world 
average Br(B ->• X sl ) = (3.52 ± 0.23 ± 0.09) • 1CT 4 [7, 8]. 

The shape of the photon energy spectrum in B — > X s j decay is sensitive to non-perturbative 
hadronic physics. At lowest order in the heavy-quark expansion, it is related to a universal 
shape function describing the momentum distribution of the b quark inside the B meson [9- 
12]. The same shape function parameterizes the leading bound-state effects in the inclusive 
semileptonic decay B — > X u l u. As a result, a precise measurement of the photon spectrum 
can be used to derive useful hadronic input for the analysis of B — > X u l v decay spectra, and 
in this way enable a precise determination of \V u b\ [13, 14]. One goal of the present paper 
is to complete the analysis of non-perturbative effects on the B — > X s 7 photon spectrum at 
subleading order in the heavy-quark expansion. This will allow us to estimate the irreducible 
theoretical uncertainties in the calculation of the B — > X s ^ branching ratio computed with 
a cut on photon energy, and it will also have implications for the extraction of \V u b\. In the 
process, we will discuss that certain terms in the standard formulae for the B — > X s 7 de- 
cay rate and photon spectrum result from an incorrect matching procedure and thus carry 
unphysical sensitivity to long-distance physics. A second goal of this paper is to properly 
factorize the short- and long-distance contributions into perturbatively calculable functions 
and non-perturbative matrix elements, using methods of effective field theory. 

The B — > Xg'-f decay rate and photon spectrum can be calculated using the optical theorem, 
which relates them to a restricted discontinuity of the forward 5-meson matrix element of the 
product of two effective weak Hamiltonians, 

dT(B -> X s7 ) oc Disc restr . [i J d 4 x (B\nl s (x)n cS (0)\B)j . (1) 

The discontinuity is restricted by the requirement that the cut propagators must include that 
of the photon and a strange quark. The effective weak Hamiltonian "H e ff consists of a sum 
of local operators, whose definitions are collected in Appendix A. The most important ones 
are the electromagnetic and chromomagnetic dipole operators Q 7j and Q 8g as well as the 
current-current operator Q\. At the lowest order in a s and l/m b , only the dipole operator 
Qj-y contributes to the decay rate. 

Note the important fact that, unlike for semileptonic inclusive S-meson decays, the B — > 
X s 7 decay rate cannot be written as the discontinuity of a forward matrix elements of time- 
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ordered products of fields. The reason is that not all cuts of the relevant Feynman graphs 
correspond to the B — > X s ~f process. For example, diagrams with penguin contractions of the 
four-quark operator Q\ contain cuts corresponding to the decay b — > ccs without a photon in 
the final state, which clearly do not contribute to the decay rate in (1). As a result, the fields 
belonging to the B — >■ X s ^ amplitude are time-ordered, while those belonging to the complex 
conjugate amplitude are anti-time-ordered. A path- integral method for the evaluation of the 
cut diagrams contributing to expressions such as (1) is the Keldysh (or time- loop) formalism 
[15, 16]. Here we will not expose the technical details of this approach (see [17] for a concise 
recent discussion), but we will mention at the appropriate places in our discussion where the 
anti-time-ordering of fields is important. 

Theoretical calculations of the forward scattering amplitude utilize the fact that Aq C d «C 
m b to express the decay rate and the photon spectrum as a series of operator matrix elements 
suppressed by powers of l/m b [18-20]. The photon spectrum has been measured accurately 
for energies E 1 > 2 GeV, and some less accurate data is available in the range between 1.7 and 
2GeV [21-25]. The partially inclusive rates obtained experimentally are defined as integrals 
over the endpoint region E < E 1 < M B /2. The shape of the photon spectrum in the region 
above 2 GeV is most useful for extracting information that can be used to determine \V ub \ 
from B X u l v decay distributions [13, 14]. Note that in the relevant region of phase space 
the variable A = m b — 2Eq is a hadronic scale of order Aqcd, which is much smaller than the 
hard scale m b of the process. In this "endpoint region" , the hadronic final state X s has large 
energy Ex ~ m b but small invariant mass Mx ~ y/m b A ~ ^m b AQ C D- This follows from the 
fact that M\ = Mb[Mb — 2-E-y), which implies that the photon energy spectrum contains the 
same information as the hadronic invariant mass distribution in B — > X s 7 decay. In this case 
the appropriate theoretical description of hadronic effects involves an expansion of the forward 
scattering amplitude in non-local operator matrix elements called shape functions [9, 10]. If 
in the future it will be possible to lower the photon cut to a value such that m b ^> A ^> Aq CD 
(this will require Eq < 1.6 GeV or so), then many (but not all) of the non-local matrix 
elements can be expanded in matrix elements of local operators using a multi-scale operator 
product expansion, consisting of a double expansion in powers of Aqcd / rn b and Aqcd / A [26] . 
However, even in the hypothetical limit E — > some non-local effects remain, which cannot 
be described using a local heavy-quark expansion in power of AQ CD /m b . 

In this paper we perform a comprehensive study of the factorization properties of the 
B — > Xg'j photon spectrum in the endpoint region, which is more systematic than previous 
analyses. Using methods of effective field theory, we propose a novel factorization formula 
valid at any order in the l/m b expansion, which is a generalization of the familiar soft-collinear 
factorization formula [27-29, 31-33] 

dT{B -+ XJP) = £ ^ £ H^4 n) ® S?° (2) 

(Tit 

n=0 b i 

for the differential distributions in the inclusive semileptonic decay B — > X u l v. Here are 
hard functions parameterizing physics at the scale m b , are jet functions describing the 
physics of the hadronic final state X u with invariant mass Mx ~ a/^Aqcd, and are soft 
functions incorporating hadronic physics associated with the scale Aqcd- The soft or shape 
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Figure 1: Graphical illustration of the three terms in the QCD factorization theorem (3) 
for B — > Xg'-f decay in the endpoint region. The dashed lines represent soft interactions, 
which must be power expanded and factored off the remaining building blocks to derive 
factorization. 



functions are defined in terms of forward matrix elements of non-local HQET operators on the 
light cone. The symbol <g> implies a convolution, which arises when the soft and jet functions 
share some common variables. 

The new element, which makes the analysis of B — > X s ^ decay more involved than that 
of semileptonic decays, is the presence of "resolved photon" contributions, which contain 
subprocesses in which the photon couples to light partons instead of connecting directly to 
the effective weak-interaction vertex [34-39]. As we will show, these subprocesses probe the 
hadronic substructure of the photon at a scale of order y / 2i? 7 AQCD- The corresponding effects 
can be described by introducing new jet functions j\ n ^ . There is no analog of this phenomenon 
in semileptonic decays, because a lepton-neutrino pair can only couple to light partons via W- 
boson exchange. The factorization formula we obtain for the photon spectrum in the endpoint 
region is 

oo 1 

dT(B -+ X sl ) = J2^H H ^jt } ® S?° (3) 

n=0 b i 
oo 1 

+ E A E H * n)j * n) ® ^ (n) ® ^ + E H i n) J i n} ® si n) ® j\ n) ® 4 n) 



n=l b 



It contains "direct photon" contributions of the same form as (2), accompanied by single and 
double resolved photon contributions that are new. Our notation is symbolic; objects denoted 
by the same symbol in the various terms refer, in general, to different quantities. Note the 
important fact that the new contributions appear first at order l/m^ in the heavy-quark 
expansion. While the jet functions j\ n ^ are cut propagator functions dressed by Wilson lines, 

~ (n) 

the jet functions J\ are given in terms of full propagator functions dressed by Wilson lines. 
A graphical illustration of the factorization formula is shown in Figure 1. 

When the photon spectrum is integrated over an interval much larger than the endpoint 
region, the direct photon contributions simplify to a series of hard coefficients multiplying 
forward £>-meson matrix elements of local operators, in analogy to what happens in semilep- 
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tonic B — > X u lP decay [18-20]. In particular, it follows that the corrections of first-order in 
Aqcd /n^b integrate to zero, since there does not exist a local, gauge-invariant operator that 
could account for such terms. An important result of our analysis is that the resolved photon 
contributions do not reduce to matrix elements of local operators in that case. Their effects 
on the total decay rate must still be described in terms of non-local operator matrix elements, 
as illustrated with a specific example in [39]. 1 

Resolved photon contributions can only arise from operators in the effective weak Hamil- 
tonian that do not contain the photon field as part of the effective, local weak interactions. 
The most important such operators are the chromomagnetic dipole operator Q 8g and the 
current-current operator Q\. Power-suppressed contributions from other operators can be 
safely neglected for phenomenological purposes. It follows that double resolved photon con- 
tributions can only arise from the operator pairs Q 8g — Q Sg , Q± — Q±, and Q\ — Q 8g , while 
single resolved photon contributions can also arise from the pairs Q$ g — Qj 1 and Q\ — Qi 1 . 
Direct photon contributions can arise from all operator pairings. Some effects involving the 
conversion of the photon into light partons have been discussed previously in the literature 
[34-39], and it is known (though not widely appreciated) that they fall outside the realm of 
the local operator product expansion. Let us comment on the various effects one by one: 

• The perturbative analysis of the Q$ g — Qsg contribution gives rise to IR-singular contri- 
butions, which at one-loop order can be regularized by introducing a non-zero mass for 
the strange quark [41]. It was argued in [34] that these singularities can be absorbed into 
the photon fragmentation functions of the strange quark and the gluon. We find that 
this factorization no longer holds in the endpoint region. Instead, the IR singularities 
must be factored into a subleading four-quark shape function. 

• The current-current operator Q\ can induce penguin-type transitions, in which two 
charm or up quarks convert into a photon and soft gluon. Previous studies of this 
effect have focused on its contribution to the total decay rate, which arises from its in- 
terference with the matrix element of Qj y [35-38]. In the present work we will generalize 
this analysis to the case of the photon spectrum in the endpoint region. 

• The square of the charm-penguin amplitude, the Q\ —Q\ double resolved photon contri- 
bution, has not yet been analyzed in the literature, but it is sometimes mentioned as a 
potentially large source of power corrections due to the fact that the operator Q\ has by 
far the largest Wilson coefficient in the effective weak Hamiltonian. We will show that 
this contribution arises first at order 1/m^ in the heavy-quark expansion. Its effects on 
the decay rate and spectrum are therefore strongly suppressed. The same is true for the 
Qi ~ Qsg double resolved photon contribution. 

• Resolved photon contributions from the Q$ g — Q77 interference term were first studied 
in [39], again with regard to their impact on the total decay rate. We will complete this 
study and generalize it to the case of the photon spectrum. 

1 The total B — > X s -f decay rate is not an infra-red (IR) safe observable. What is usually meant by this 
term is the rate defined with a very low cut on photon energy, and with a subtraction of duality-violating 
charmonium resonance contributions [40]. 
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We begin our analysis with a review of known results for the B — > X s ^ photon spectrum 
in Section 2, indicating a couple of problematic features in the formulae routinely used in 
the literature. The new factorization formula (3) will be derived in Section 3 using a two- 
step matching procedure from QCD to soft-collinear effective theory (SCET) [28, 42, 43] 
and heavy-quark effective theory (HQET) [44]. In Section 4 we discuss the factorization 
properties of the various contributions to the B — > X s ^ photon spectrum, which arise from 
different pairs of operators in the effective weak Hamiltonian. This includes, in particular, a 
detailed discussion of the new subleading shape functions for the contributions from operator 
pairs other than Q 77 — Q 7l . These have not been considered previously in the literature, 
except for a particular subleading shape-function contribution to the total B — > X s ^ decay 
rate arising from the operator pair Q 7l — Q$ g [39]. In Section 5 we use the invariance of 
the strong interaction under the discrete symmetry PT to prove that the subleading soft 
functions are real, i.e., they do not carry non-trivial strong phases. The implications of our 
findings for the integrated B — > X s ^ decay rate are studied in Section 6. We show that the 
resolved photon contributions must still be described in terms of matrix elements of non- 
local operators, whose effects cannot be reduced by lowering the cutoff on the photon energy. 
Finally, in Section 7 we study the phenomenological implications of our results by estimating 
the irreducible theoretical uncertainty in the prediction for the B — > X s 7 decay rate integrated 
over the range E 1 > 1.6 GeV. We then summarize our results and give some conclusions. 
Three Appendices contain our definitions of the operators in the effective weak Hamiltonian, 
a summary of input parameters, and a detailed exposition of the matching of the effective 
weak Hamiltonian onto operators in SCET. Readers not interested in the technical details of 
our derivations should consult Section 2 and then proceed with Sections 6 and 7. 

Even though we only sketch the derivation of the new factorization formula in Section 3, 
we consider this discussion as solid as that for many other processes discussed in the context of 
SCET. Still, we do not claim to have a rigorous proof of factorization. Indeed, in our analysis 
we will encounter one particular contribution to the B — > X s 7 photon spectrum and decay rate, 
for which the resulting convolution integrals derived using SCET suffer from an ultra-violet 
(UV) divergence. In its current formulation, the effective theory does not provide a systematic 
framework for regularizing this divergence. The problem of divergent convolution integrals in 
SCET has been encountered previously in the context of heavy-to-light form factors [45-48] 
and power-suppressed contributions to hadronic fi-meson decays [49, 50]. It is to some extent 
still an open question whether these integrals indicate a failure of factorization, or whether 
they can be cured by a generalization of the theoretical framework of SCET (an attempt in 
this direction was initiated in [51]). An important difference is that in all previous cases these 
divergences were of IR origin. In our case, the convolution integrals diverge in the UV. Such 
divergences appear to be rather generic in the description of higher-order power corrections, 
because the resulting convolution integrals contain higher powers of soft momentum variables. 
The physical origin of the divergence and its interpretation are entirely transparent. Still, 
the presence of this effect is problematic for the consistency of SCET as a bona fide effective 
field theory and calls for a cure. We will discuss a simple treatment of the divergence using 
a hard cutoff on the convolution integrals. We do not claim, however, to have a systematic 
procedure that would work at higher orders in perturbation theory and allow for a consistent 
resummation of large logarithms. In that sense our derivation of factorization is incomplete. 
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2 Review of known results and preview of new ones 



Let us briefly summarize what is known in the literature about the various terms in the 
factorization formula (3). Separating the contributions from different operators in the effective 
weak Hamiltonian, we write the heavy-quark expansion of the CP-averaged B — > X s ^ photon- 
energy spectrum in the endpoint region p + = mj — 1E 1 = C(Aqcd) in the form 



dT G%a\V t ,V*\ 2 



F tX\V tb V ts \ 2 3 



dE 1 2tt 4 



\H y (fi)\ 2 / dum b j(m b (u + p + ),fi)S( y u,fi) 
J-P+ 

( 4 ) 

+ — J2 Re [ C *M CM] M&v v) + ■ ■ 



where A = MB—m b , and the ellipses represent terms of order 1/m 2 and higher. For convenience 
we have factored out two powers of the running 6-quark mass (defined in the MS scheme) and 
three powers of the photon energy, as this is the correct energy dependence of the leading 
contribution to the spectrum. 

The term in the first line of (4) is the leading-power contribution and is well understood 
theoretically. At this order the effective weak Hamiltonian for B — > X s 7 decay matches onto 
a unique leading-order current operator in SCET. The hard matching coefficient H^(fi) = 
CV 7 (yu) + 0(a s ) for this current receives contributions from all operators in the effective weak 
Hamiltonian, not just Q7 7 , as soon as one goes beyond the leading order in perturbation theory 
[26]. The contribution proportional to CV 7 is known to order a 2 [52, 53], while the remaining 
terms are known to order a s . When the effective current operator is further matched onto 
HQET, a single jet function J(p 2 , /i) = 5{p 2 ) + (D(a s ) arises, which is given by the discontinuity 
of the quark propagator in light-cone gauge and has been calculated to two- loop order in [54] . 
The remaining HQET matrix element defines a single, leading-order shape function via [10] 

S(uj U )= f- e-** (B(v)\h(tn)S n (tn)St(0)HO)\B(v)} 
{ ' W J 2tt 2M b ' 

Here v denotes the four-velocity of the B meson, and n is a light-like vector pointing along the 
direction of the final-state hadronic jet. We normalize these vectors such that v 2 = 1, n 2 = 0, 
v ■ n = 1, and v° > 1. The soft Wilson line S n is defined as 



S n (x) = P exp f ig j dun ■ A s (x + un)\ , (6) 



where the path-ordering symbol P means that fields with larger u values stand to the left of 
those with smaller ones. The conjugate Wilson line has the opposite ordering prescription. 
These definitions imply that S' n (tn)S'^(0) = [tn, 0] is a straight line segment connecting the 
points tn and 0, with gauge fields closer to the point tn standing to the left of those closer to 0. 
Taking the complex conjugate of relation (5) and using translational invariance, one finds that 
the shape function is real. The functions H 7 , J, and S incorporate contributions associated 
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with different scales in the problem. The hard function if 7 receives virtual corrections of 
order the hard scale fih ~ n^b, while the shape function S encodes non-perturbative hadronic 
physics associated with the soft scale fi s ~ p + ~ Aqcd- The jet function describes the 
properties of the final-state hadronic jet, whose invariant mass scales like fihc ~ a/ Aqcd in 
the endpoint region. This intermediate scale is the scale of (anti-)hard-collinear virtualities. 
Large logarithms arising from ratios of these various scales can be resummed to all orders in 
perturbation theory by solving renormalization-group equations in the effective theory [26-28] . 

Beyond the leading power in the heavy-quark expansion, the proper factorization of the 
various contributions to the decay rate has not yet been discussed systematically in the liter- 
ature. Such an analysis is the main goal of the present work. In phenomenological discussions 
of B — )> X s 7 decay one usually starts from expressions for the power-suppressed terms de- 
rived in the naive parton model, i.e., by computing the inclusive decay of an on-shell b quark 
[12, 41, 55]. Including only the phenomenologically relevant contributions from operator prod- 
ucts of Ql, Q-j-y, and Q$ g , and setting V u b = for simplicity, this yields for the first-order power 
corrections 

Fr (^) = (l6 In ^ - ISV 

pparwp x _ C F a a (n) (2 m bP+ 1 



CfoJjj) 10 / 7 \ 
4tt 3 ' 



.part/p C F a s (fi) f 2\ m b \ 2 



In this paper we adopt the scaling m 2 c = 0(rribAQCD) for the charm-quark mass, meaning that 
the ratio m 2 c /mb remains a constant of order Aqcd in the heavy-quark limit. Note that the 
expression for Ffj T includes a non-perturbative effect proportional to the HQET parameter 
A 2 = — M|)/4 0.12 GeV 2 [35-38], which is of the same order in power counting as the 
perturbative contribution. It is related to charm-penguin diagrams with a soft gluon emission. 
If we were to adopt the alternative counting scheme where m c = 0(mb) in the heavy-quark 
limit, then some of the expressions in (7) would change. In that case 



*n^)-^jf*(i->|i-^) 

Fn^c) = -3i^(S 7 ,„) = f-i) fjxxRe [l - F(j) 

where z = (nic/mb) 2 , and we have defined the penguin function 



(8) 



F(x) = Ax arctan 2 ( ) . (9) 
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The non-perturbative contribution to F^ rt would be power-suppressed in that case and should 
be dropped for consistency. 

In order to account for non-perturbative effects other than those described by the A2 term, 
the simplest recipe used in the literature is to replace p + — > u> + p + in the above expressions 
and convolute them with the leading-order shape function S(u,fi), e.g. 

F 77 ( Ej , $ = £ du (l6 In - is) n) , 

u P+ (10) 

and similarly for the other terms [12, 13]. 

Beyond the leading order in l/m b , the photon spectrum also receives contributions in- 
volving more complicated soft functions, usually called subleading shape functions. So far 
they have been studied only for the direct photon contribution from two insertions of the 
electromagnetic dipole operator Q 7l [30-33]. In the notation of [32], one obtains at tree level 

F 7 s 7 SF (£ 7 ,/i) =p + S(-p + ,fi) + s(-p+,fi) -t(-p + ,fi) +u(-p+,fi) -v(-p + ,fi) 

[fi s \-P+^) + fi s) (-P+,ri] 

These same functions also contribute, in other combinations, to the semileptonic B — > X u l v 
decay spectra in the endpoint region. One can think of t and v as non-local generalizations 
of the _B-meson matrix element of the subleading HQET chromomagnetic operator, and of 
11 as a non-local generalization of the matrix element of the kinetic operator. The function 
s arises from an insertion of the subleading HQET Lagrangian into the matrix element for 
the leading shape function in (5). The functions fffi and /® arise from the matrix elements 
of non-local four-quark operators. In the present paper we will encounter other four-quark 
shape functions, which are unique to radiative decays. It is therefore hopeless to try to find 
weighted distributions of B — > X s 7 and B — > X u l v events, in which the subleading shape 
functions enter in the same combinations - a goal pursued in [56], working at tree level and 
neglecting all operators in the effective weak Hamiltonian except Q 7l . As in all previous anal- 
yses of subleading shape-function contributions, it is sufficient for phenomenological purposes 
to restrict the analysis to the tree level, since so little is known about the functional forms of 
the subleading shape functions. In the language of the factorization formula (3), this means 
that the corresponding hard and jet functions are computed at zeroth order in a s /ir. We 
do, however, include jet functions associated with a factor g 2 = 4ira s , which can arise from 
tree-level hard-collinear gluon exchange. In this work we complete the analysis of subleading 
shape functions for B — > X s ^ decay by analyzing the contributions analogous to (11) for the 
remaining pairs of operators in the effective weak Hamiltonian. 

It would be incorrect to simply add the partonic contributions and the contributions from 
subleading shape functions, such as (10) and (11), as this would lead to double counting. In 
fact, since the partonic expressions (7) and (8) have not been derived from a proper matching 
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procedure, they secretly contain some soft contributions, which should be subtracted and 
absorbed into the subleading shape functions. The explicit expressions for F^ rt and F^ r 
in (7), which contain parametrically large logarithms, already hint at the fact that such a 
subtraction is required. The dependence of Fg 8 ar on the strange-quark mass is clearly a sign 
of an unphysical sensitivity to the IR region, which should not be present in a short- distance 
coefficient function. For the case of F^ Tt one might think that the large logarithm results from 
a combination of hard-collinear and hard scales, ln(mft/p+) = ln(m^) — ln(m&p+), in which case 
it would have a short-distance origin. We will see, however, that it results from a combination 
of hard-collinear and soft scales, \n(m b /p + ) = \n(rribP+) — 2 ln(p + ). The sensitivity to the soft 
scale p + must be subtracted and absorbed into a subleading shape function. 
Our improved expressions for the coefficient functions read 

F 77 (E„ fx) = £ du (l6 In mb( ^+ P+) + 9 ) S(u>, fx) + iff \E„ fx) , 

F 88 (£ 7 ,/i) = CfC ^-^ j duj Q ln m b(u+P+) _ 1^ $(uj, fx) + 47ia s (fx) f 88 (-p + , fx) , 
F 78 (E 1 ,fx) = ^^-f £ cLoS(u : fx) + 4na s (fx)Re[f$(-p +: fx) + fW (12) 
F 17 (E 1:t x) = f-^j f dwS(u : fx)+J2 S q Ref 17 , q (-p + ,fx), 



q=c,u 



F u (E,,fx) = F 18 (E,,fx) = ^^l j dcoS^fx), 



p+ 



5 « = 1 l . x « = w . ( i3 ) 



where we now also include the effects of up-quark loops in Fi 7 , i.e., we no longer set V u b 
We have defined 

|A,| 2 Re [Ci(a*)C?»] 

where 5 C + 5 U = 1 by unitarity of the CKM matrix. Note that 5 U is of second order in the 
Wolfenstein parameter A ~ 0.22 and thus numerically very small. The perturbative terms 
involving convolutions with S(u, fx) are now free of IR-sensitive contributions, and all long- 
distance physics resides in the shape functions. 

We finish this section with an important remark, which is somewhat orthogonal to the 
main thrust of this paper but nevertheless relevant. In the analysis of the B — > X s 7 decay 
rate and photon spectrum, it is customary to adopt for the charm-quark mass a running mass 
defined at a hard-collinear scale jihc ~ a/"^Aqcd ~ m c [57]. For instance, the default value 
adopted in [5] is m c = m c (1.5GeV). This scale choice is indeed appropriate for charm-quark 
mass effects residing in the jet functions entering the factorization formula (3). We will be 
concerned with such effects in Section 4.3. On the other hand, charm-quark mass effects also 
enter some of the hard functions in the factorization formula, for instance via the coefficient 
Hy(fx) in (4) [26], or via phase-space functions such as those shown in (8). In this case the 
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charm-penguin loops are probed at virtualities of order m b , and it is therefore appropriate to 
use a running mass m c = rn c (nh) evaluated at a hard scale fih ~ n^b- This can have important 
numerical effects, enhancing the theoretical prediction for the total decay rate by up to 3%. 

3 Schematic derivation of the factorization formula 

The endpoint region of the B — > X s ^ photon spectrum is defined as the kinematical region 
where the hadronic jet X s has large energy compared to its invariant mass: E x ~ rn b , but 
M x ~ y^m^AQCD- We define two light-like vectors, n and n, which are aligned with the 
directions of the hadronic jet with total momentum P x and the photon with momentum q. 
These two vectors satisfy n + n = 2v, where v is the velocity of the B meson. Specifically, 
we have P x = E x n + (M X /AE X ) n and q = E^n. In the 5-meson rest frame, we may choose 
n M = (1, 0, 0, 1) and = (1, 0, 0, —1). It is convenient to decompose 4- vectors in the light-cone 
basis spanned by n and n, 

a M = n-a—+n-a— + a^ = a% + a^_ + a^. (14) 

We will often use the short-hand notation a ~ (n ■ a,n • a,aj_). Note that by definition the 
external momenta P x , q, and MbV have vanishing perpendicular components. 

SCET and HQET are the appropriate effective field theories to study the factorization 
properties of inclusive I?-meson decay spectra in the endpoint region [28, 29, 31-33]. We will 
need several types of SCET modes for our analysis, each one corresponding to a particular 
physical scale relevant to the process. Indeed, once the relevant modes have been identified, 
factorization follows from a sequence of simple and by now familiar steps. High-energy scales 
such as the electro- weak scales m t , My/ or the hard scale set by the heavy-quark mass m?, 
are integrated out before one enters the low-energy effective theories, and hence there are no 
fields for such hard quantum fluctuations in SCET or HQET. 

The expansion parameter in the factorization analysis is A = AQ CD /m 6 , where we do 
not distinguish between M#, nib, E x , and E 7 , all of which are of the same order in the 
endpoint region. The partons that make up the final-state hadronic jet X s carry momenta 
that generically scale like the total jet momentum P x . In light-cone components this implies 
Phc ~ m b (A, 1, \/A). We will refer to modes with this scaling behavior as hard-collinear (he) 
fields. The final state as well as the initial B meson also contain soft partons with momenta 
Ps ~ m b (A, A, A) of order Aq CD . The light-like photon momentum q itself does not set a 
physical scale; however, partons with momenta scaling like pj^ ~ m fe (l,A,vT) can convert 
into a photon accompanied by soft partons, which can be absorbed by the hadronic final state 
or originate from the initial B meson. We will refer to modes with this scaling behavior as 
anti-hard-collinear (he). The invariant mass of a set of anti-hard-collinear partons scales like 
y^m^AQCD- Note that processes in which a parton fragments into a photon plus an energetic 
parton moving along the n direction, such as those considered in [34], are kinematically not 
allowed in the endpoint region. The produced energetic parton cannot be absorbed by the 
low-mass final-state hadronic jet. 
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In this paper, we denote by ^ c and £^ the effective SCET fields for hard-collinear and 
anti-hard-collinear quarks, respectively. We define them as 

i h c = wH n ~V\, fc = wte ii ~>/\, (15) 

where £ n and £ s are two-component spinor fields in the effective Lagrangian, which obey 
= f££n = 0. The quantities W n and W n are the familiar (anti-)hard-collinear Wilson lines 
of SCET. Similarly, we define the hard-collinear and anti-hard-collinear gluon fields [58] 

Kc = WJ (iKW*) ~ ( A > °> v^A) , A^= W\ (iD^W n ) ~ (0, A, y/X) . (16) 

Finally, we need the soft heavy- and light-quark fields h,q ~ A 3 / 2 and the soft gluon field 
~ (A, A, A). The effective heavy-quark field obeys jh = h. For clarity, soft light-quark fields 
will often be denoted by the flavor of the corresponding particles (q = u,d,s, . . .). Finally, 
note that in general the gluon fields have the same scaling properties as the corresponding 
momenta (apart from the large components of the (anti-)hard-collinear gluon fields, which 
have been "gauged away" by the introduction of the Wilson lines), and the same is true for 
derivatives acting on these fields. 

The effective fields defined in (15) and (16) are invariant under (anti-)hard-collinear gauge 
transformations, while they transform homogeneously under soft gauge transformations [43]. 
An important technical step in deriving factorization formulae using SCET is the decoupling 
of the soft gluons from the (anti-) hard-collinear fields. This is accomplished by introducing 
the soft Wilson lines S n (x) in (6) and corresponding Wilson lines S n defined analogously with 
n replaced by n. We then perform the decoupling transformations [28, 59] 

&c(s) = S n (x.) , Al(x) = S n (x_) A^(x) Sj(x_) , (17) 

and similarly for the anti-hard-collinear fields. When expressed in terms of the "sterile" fields 
with superscript "(0)", the SCET Lagrangian no longer contains interactions between soft and 
(anti-)hard-collinear fields at leading order in A. 

The derivation of the factorization formula (3) proceeds as follows. In the first step, the 
QCD matrix element of the product of the two effective weak Hamiltonians in (1) is matched 
onto operators in the effective theory SCET(/ic, he, s), which is the version of SCET containing 
hard-collinear, anti-hard-collinear, and soft degrees of freedom. In this step the hard functions 
appear as Wilson coefficients capturing the effects of hard quantum fluctuations. A priori, 
an operator in the effective weak Hamiltonian can be matched onto any operator in SCET 
with the right quantum numbers. One finds, however, that up to order l/m& only a small 
number of effective-theory operators contribute. Many aspects of this first matching step have 
been discussed in [60], where a factorization theorem was derived for the exclusive decays 
B — > fT*7 and B — > p^y. Note that in the endpoint region anti-hard-collinear partons cannot 
be part of the final state X s , because this would lead to an invariant hadronic mass Mx ~ ra^, 
in contrast with the required scaling Mx ~ a/^&Aqcd- We therefore need insertions from the 
SCET Lagrangian, which convert all anti-hard-collinear particles into a photon plus a set of 
soft partons. As will be discussed in Section 4, these insertions are always power suppressed. 
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After this is done, any given contribution to the decay rate is related, in position space, to a 
matrix element of the generic form 



Note that we have introduced two types of anti-hard-collinear fields. The fields 0^ (<f>j-) only 
couple to the photon connected to the initial (final) 5-meson state via the weak effective 
Hamiltonian, see Figure 1. These fields should be thought of as different entities, which do 
not interact with one another (i.e., there are no interaction terms in the effective Lagrangian 
coupling (fij^; and 0^)- The reason is that the exchange of an anti-hard-collinear particle 
across the two sides of the forward amplitude is forbidden kinematically, as this would lead to 
a hadronic final state with mass of order m^. In other words, anti-hard-collinear propagators 
are never cut. 

The soft heavy-quark fields need to be part of any effective-theory operator. The three 
brackets [. . . ] can contain generic products of soft, hard-collinear, and anti-hard-collinear fields. 
The presence of a hard-collinear jet in the final state requires that the hard-collinear bracket 
must contain at least two fields, either a pair of strange quarks or of other light partons (gluons 
or quarks). In the latter case, the strange-quark pair must appear in the soft bracket. The 
anti-hard-collinear bracket, on the other hand, can be empty. Recall that any field in the 
effective Lagrangian scales like a positive power of V% so adding an additional field to an 
operator always leads to further power suppression. 

Because of the particular scaling of soft and (anti-)hard-collinear momenta, the soft fields 
must be multipole expanded when they couple to (anti-)hard-collinear fields [61]. The mul- 
tipole expansion is subtle if there are external momenta in the problem whose components 
are nearly identical, such as n ■ (rribv) n • q in our case. The correct form of the expansion 
must then be determined on a case-by-case basis for each operator, rather than derived from 
a simple set of rules. We find that the heavy-quark fields must always be expanded about x_ , 
while other soft fields can depend on either x_ or x + , or both. 

At this point the fields belonging to the three types of modes can still interact with each 
other by the exchange of soft gluons. These interactions are unsuppressed in SCET, but they 
have an eikonal structure and can be removed by field redefinitions [28, 59]. The remaining 
power-suppressed interactions are treated as Lagrangian insertions and so are included as parts 
of the operators in (18). After the decoupling transformations, the forward 5-meson matrix 
elements needed for the calculation of inclusive decay spectra can therefore be factorized into a 
S-meson matrix element of soft fields multiplying vacuum expectation values of hard-collinear 
and anti-hard-collinear fields: 



In this process the soft Wilson lines from (6) arise, which must be included as part of the soft 
matrix elements. In fact, they render these non-local matrix element gauge invariant. 

In the last step, we match SCET onto HQET by integrating out the (anti-)hard-collinear 
fields. This can be done perturbatively, because the corresponding scales are in the short- 
distance regime. The Wilson coefficients of this matching are simply the perturbative expres- 
sions for the vacuum correlation functions of the (anti-)hard-collinear fields. Their Fourier 



Disc {B(v)\h(x)[Uxi) ■ ■ ■ ]h(0) frUVj) ■■■} i<k-c(z k ) ■ ■ ■ <j>' r M) ■ ■ -Wv)) . 



(18) 



(B(v)\h(x-)[M*i*) ■ ■■M0)\B(v)) x Disc (Op^) . . .]|0> 
x(O|[0£W..]|O> (O|0g>(*0...]|O>. 



(19) 
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transforms define the momentum-space jet functions J- n ^ and J\ n> , and the Fourier trans- 



forms of the soft matrix elements define the soft functions S} 



J 



(n) 



Disc(0|[tf c ) ( % )...]|0) 



(n) 



F.T. 



<0|[^(z fc )---]|0> 



F.T. 



7 (n) 



{B\h(x^)[<p s (x iT ) . ..}h(0)\B) 



(20) 



J F.T. 



If both anti-hard-collinear brackets in (19) are empty, the corresponding contribution to the 
spectrum becomes part of the first term in the factorization formula (3). If only one of them is 
empty, the contribution becomes part of the second term. Finally, if both anti-hard-collinear 
brackets are not empty, the contribution belongs to the last term. Note that the momentum- 
space functions defined above will, in general, depend on several Fourier variables (one less 
than the number of space-time variables in the corresponding non-local matrix elements), and 
as a result the convolutions in the factorization formula (3) involve multi-dimensional integrals. 

This concludes the derivation of the factorization formula. Our main focus in this work is 
on the resolved photon contributions giving rise to the second and third terms in (3). We will 
now present a detailed discussion of the matching procedure for these contributions. 



4 Systematic analysis of resolved photon contributions 

4.1 Matching onto SCET 

As outlined in the previous section, the analysis of the B — > X s ^y photon spectrum in the 
effective theory consists of two steps. In the first step the effective weak Hamiltonian sum- 
marized in Appendix A is matched onto operators in SCET consisting of (anti-)hard-collinear 
and soft fields. The Wilson coefficients arising in this step are the hard functions in 
(3). In the second step the (anti-)hard-collinear modes are integrated out, and the theory is 
matched onto HQET. The Wilson coefficients arising in this step are the jet functions and 
j\ n) . The remaining hadronic matrix elements define the soft functions S^ n \ For simplicity, 
we will restrict ourselves to the tree-level approximation for hard quantum corrections, and to 
the one- loop approximation for (anti-)hard-collinear quantum fluctuations associated with the 
leading shape function in (5). The Wilson coefficients of the new subleading shape functions 
will be computed at tree level, but including contributions of order g 2 = Ana s resulting from 
tree-level (anti-)hard-collinear gluon exchange. 

In constructing the possible operator basis of SCET for B — > X s ^ decay, we require that 
the final state should contain only one anti-hard-collinear particle, which is the photon field 
with momentum = E^n. All the other particles in the final state, including one strange 
quark, need to be either hard-collinear or soft. At least one of the particles in the final 
state must be hard-collinear, either the strange quark or a gluon. Note that we can have 
several anti-hard-collinear and/or soft fields be part of the possible operators, provided that 
all the anti-hard-collinear particles are converted into the photon plus soft particles via SCET 
Lagrangian insertions. The number of soft particles in the final state is restricted only in 
the sense that adding soft fields to an operator always leads to power suppression. From 
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Figure 2: 0(A 1/2 ) (top row) and (9(A) (bottom row) conversions of anti-hard-collinear 
particles into a photon accompanied by soft particles. Only some representative dia- 
grams are shown. 



these simple requirements, it is straightforward to find the possible SCET operator basis 
systematically. 

Since the SCET expansion parameter is v^A with A ~ Aqcd/?™&, we need to consider oper- 
ators of leading, next-to-leading, and next-to-next-to-leading order in SCET power counting 
in order to systematically analyze Aqco/^ft-suppressed contributions to the photon spectrum 
[43, 47, 60]. The possible effective weak-interaction operators can then be divided into two 
classes, depending on whether they contain a photon field or not. The first class of operators 
leads to direct photon contributions, whereas the second class gives rise to resolved photon con- 
tributions. In the latter case, the operators must be supplemented with Lagrangian insertions 
that convert the anti-hard-collinear particles into a photon. 

As our main interest is in the resolved photon contributions, we begin by systematically 
analyzing how anti-hard-collinear partons can be converted into a photon. These conversions 
can be derived from the SCET(/ic, he, s) Lagrangian. For power counting the photon field 
scales like an anti-hard-collinear field, and hence the conversion of any number of anti-hard- 
collinear particles into a photon is unsuppressed as long as it is allowed by the rules of the 
leading-order SCET Lagrangian. However, each conversion involving a soft parton costs a 
certain power of \/X. At 0{\f\) we find the possibilities (here and below one could add any 
number of anti-hard-collinear gluons on the left-hand side of the relations) 

fa^A? + q, fc + fc^A? + A., Atz + Ak^AT + A.. (21) 

Similarly, at 0(A) we have 

■%-^r + 9 + 9, A 1 z^A? + A. + A s . (22) 

These five possibilities are illustrated in Figure 2. In each case the last conversion, which 
involves three gluon fields, is not needed for our tree-level analysis. 

In the matching of the effective weak Hamiltonian onto SCET we focus on the contributions 
of the operators Q77, Qsg, and Q c { u . Other four-quark operators, which could be treated in 
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a way analogous to Ql' u , give rise to negligible effects. There is a large number of SCET 
operators appearing in the matching relations, all of which are listed in Appendix C. Here 
we will explicitly present only those operators that are needed for our tree-level analysis of 
resolved photon contributions. The unique leading-order operator arising in the matching 
relation for Q 7l is 

Q 7l (x) ^ e**™ l hc {x) | [in ■ d4?(x)] (1 + 75) h(x.) , (23) 

where the derivative acting on the photon field produces a factor —n-q = —2E^. This operator 
arises at 0(\ 5 ^ 2 ) in SCET power counting. It is shown in the first row in Figure 3. The power 
suppression of other operators must be evaluated in comparison with this scaling. While at 
tree level only Q 7l matches onto the operator in (23), at 0(a s ) the other operators contribute 
as well, with Wilson coefficients that can be found in [26]. From the many operators arising at 
subleading power in SCET, we only need those that can give rise to l/m^-suppressed resolved 
photon contributions. They must contain at least one anti-hard-collinear field. From the 
matching relation for Q 8g , we need the leading-order operators 



ihc(x) - [in ■ d4^ c± (x)} (I + 75) h(x-) 



+ feO*0 £ [in ■ dJfi hcl _{x)\ (1 + 75) h{x.) 



(24) 



They are shown in the second row in Figure 3. The conversions of the anti-hard-collinear fields 
into the photon plus soft fields follow from (21) and (22) and give rise to power suppression, as 
indicated in the figure. Finally, from the matching relation for the current-current operators 
Ql, we need the C(A 3 ) operator 

Ql(x) fcfil - lb )h{x-) fe(x) 7M (l - 75)^) (25) 

illustrated in the third row in Figure 3. Here ^ c is the strange-quark field, while ^ are fields 
for the quark with flavor q = c, u in the four-quark operator. The anti-hard-collinear quark 
pair must be converted into a photon plus a hard-collinear or soft gluon via a penguin loop, 
as indicated by the second relation in (21) and the second graph in the first row of Figure 2. 

In our analysis we will also include power corrections of 0(a s ) to the direct photon con- 
tributions in the factorization formula (3). They involve a hard-collinear loop and give rise 
to a convolution of subleading jet functions with the leading shape function. The operators 
Qs g , and Ql' u can all be matched onto 0(\ 3 ) and 0(\ 7 ^ 2 ) SCET operators containing 
a photon field accompanied by hard-collinear quark and gluon fields. A complete list of the 
corresponding operators can be found in Appendix C. A specific example, which we will need 
in the analysis in Section 4.5, is the unique 0(X 3 ) four-particle operator containing the photon 
field in the matching relation for the dipole operator Q 8g , which reads 

Q 8g (x) =^ e- im >™ i hc {x) e d eft?^) [™ • d4 hcl _{x)\ (1 + 75) h{x.) , (26) 

4?r in - d 
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Figure 3: Relevant operators arising in the matching of the effective weak Hamiltonian 
onto SCET. While many other operators exist, only those shown here contribute to 
the resolved photon contributions at tree-level in perturbative matching. 

where = —1/3 is the electric charge of a down-type quark in units of e. We refrain from 
listing the relevant operators descending from Q 7l and Q c { u here. The resulting OiXjm^) 
direct photon contributions are of the form shown in Figure 4, where the two diagrams on the 
left show the possible operator products in SCET. In the second graph, the operator on the 
right vertex represents the leading-order matching contribution given in (23). This diagram 
therefore only arises in pairings of the form Qi — Q^. The graph on the right illustrates the 
structure of soft fields remaining after the hard-collinear fields have been integrated out in the 
second matching step. Here the dashed horizontal line represents a Wilson line along the n 
direction. In this case the soft matrix element in HQET is the leading shape function S(u>, jj). 




Figure 4: Diagrams representing the OiXjm^) direct photon contributions arising from 
hard-collinear loops. The two graphs on the left represent products of SCET operators 
(see Appendix C for a complete list), while the graph on the right represents non-local 
operators built out of the soft fields remaining after the matching onto HQET. 
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4.2 Analysis of the Q 7l — Q 7l contribution 

We have discussed in Section 2 that at order l/m& in the heavy-quark expansion the Qj^ — Qi-y 
terms in the B — > X s j photon spectrum contain a parton-model contribution of order a s as 
well as a tree-level contribution from subleading shape functions. They are given in (7) and 
(11), respectively We have argued that F 7 p 7 a ' Tt (E 1 , fx) contains a soft contribution, which needs 
to be extracted and absorbed into the subleading shape-function contribution to avoid double 
counting. We will now discuss this subtraction in more detail. 

Evaluating the diagrams in Figure 4 with two insertions of Q 7l , we find that they give rise 
to a convolution of the leading shape function in (5) with a subleading jet function J su b\{p 2 ) 
resulting from the cuts of hard-collinear loop diagrams. This jet function is divergent, and 
using dimensional regularization with d = 4 — 2e space-time dimensions we obtain the bare 
expression [62] 

^ (P 2 ) = Hf) ^ ^ ("f - 16 In £ + 9 ) , (27) 

where 1/e = 1/e — 7^ + ln47r, and the factor of (d — 2) arises from the Dirac algebra in d 
dimensions. Subtracting the pole term, we obtain 

F£\E y ,ri = °*^ £ ^^i 61n ^^+^) +25 -16c R s) S( U ,ri, (28) 

where the scheme-dependent constant crs vanishes in the MS scheme, = 0. The \n(m b ) 
term in the above result agrees with that in the parton-model inspired expression for F 77 given 
in (10). However, the p + -dependent terms and the constant accompanying the logarithm are 
different in the two expressions. This difference is accounted for by the subleading shape- 
function contribution F$(E 1 ,fi) = Ff 7 SF (E 1 , /x) given in (11). 

In order to show that this is indeed the case, we derive the partonic expressions resulting 
from a naive, perturbative calculation of the subleading shape functions. To this end, we need 
to calculate the one-loop corrections to the matrix elements of the soft operators corresponding 
to the subleading shape functions uS(u,fi), s(w,/x), t(u,fi), u(uu,fi), v(uj,fj), as well as the 
tree-level matrix elements of the soft operators corresponding to f u (oj, fi) and f v (oo, fi) between 
on-shell heavy-quark states with velocity v. Without loss of generality we take v± = 0, 
and as a result the matrix elements of the soft operators corresponding to t and v vanish 
at one-loop order, since they only contain gluons with transverse polarization. The tree- 
level matrix elements of the soft operators corresponding to f u and f v also vanish, since 
they contain scaleless integrals over the n-components of the light-quark momenta. Finally, 
the matrix element of the soft operator corresponding to s is non-zero only when the heavy 
quarks are off shell with a non-zero n-component of the residual momentum. As a result it 
does not contribute when we set the residual momentum to zero in the partonic calculation. 
We are therefore left with only u S and u. The former can be extracted from the calculations 
performed in [29]. After setting v ■ k to zero, we find 

oj S barc (u,) = 9(-u) (-1 - 4 In ^ + 4) . (29) 
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For the latter, one finds by explicit calculation that [62] 



Mbarc (c) = 9{-u) + 12 In ^ - 20) . (30) 

Subtracting the UV poles in the MS scheme, we obtain from (11) 

^f F (^)L=^(«lng-24), (31) 

where the subscript "pert" indicates that these are naive perturbative expressions for non- 
perturbative hadronic functions. Adding to this result the naive perturbative expression for 
Fjf(E 1 ,ii) in (28) obtained by replacing S(u,fi) — > 5(u) yields 

^>(^)U +f?m^)U = ( 161 v + 1 - 16c -) • < 32 > 

which coincides with expression for Fj 7 Tt (E 1 ,fi) in (7) apart from the constant term. 

The remaining difference has its origin in the (d — 2) prefactor in (27), which results from 
the d- dimensional Dirac algebra. The reason is that relation (11) was derived in [32] by 
working in d = 4 dimensions. Indeed, it is convenient to renormalize the subleading shape 
functions before evaluating the traces arising from the Dirac structures specific to a given 
process. Only in that case the definitions of the subleading shape functions are the same for 
different processes, such as B — > X s 7 and B — > X u l v decays. However, in this case it is crucial 
that the same subtraction scheme is adopted in the calculation of the subleading jet-function 
contribution in (28). In other words, we should work in the DR subtraction scheme [63, 64], 
where the Dirac algebra is performed in d — 4 dimensions, while loop integrals are evaluated 
with d = 4 — 2e. In this scheme crs = 1 in (28) and (32), and the latter expression then agrees 
with that in (7). Our final result is therefore 

F 77 (E 17 „) = £ du (l6 In mfe( ^+ P+) + 9) S(u, „) + Ff F (£ 7 , „). (33) 

This is the first equation in (12), which replaces the "incorrect" (in the sense of improper 
factorization) result in (7). 

Throughout this paper we set the strange-quark mass to zero, which turns out to be an ex- 
cellent approximation numerically. Let us nevertheless briefly comment on finite strange-quark 
mass effects. Taking m s to be non-zero gives rise to a subleading jet function proportional 
to rn 2 s /p\ c) where ph c is the momentum of the hard-collinear jet containing the strange quark 
[65]. If we adopt the scaling m s ~ (9(Aq CD ), this function scales as AQ CD /m b in the endpoint 
region and contributes to F 77 . From [65], we find that the extra contribution is 

F^ a (E 1 ,fi) = -H 77 dum b j m (u + p + ,fi) S(u,fi) , (34) 
J-P+ 
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Figure 5: Diagrams arising from the matching of the Q\ — Q 7l contribution onto 
SCET (left) and HQET (right). Horizontal (vertical) dashed lines denote non- localities 
obtained after (anti-)hard-collinear fields have been integrated out. 



where we have defined j m = \ Im J™ . Both H 77 = 1 + 0(a s ) and J™ fc can be found in [65]. 
At the lowest order in a s , 

2 

in 

j m {u + p+, fi) = — S'(u + p+) + 0(a s ) , (35) 

771;, 

and indeed j m is suppressed by m 2 /(mf,AQCD) compared to the leading-order jet function 
J(p 2 ,fi). One could argue that for m s w 100 MeV and Aqcd ~ 500 MeV the parameter m s 
should scale as a higher power of the SCET expansion parameter A, e.g. m s ~ A 2 . This would 
imply that j m can be neglected at order AQ CD /m b . Contributions of j m to coefficient functions 
Fij other than F 77 are further suppressed by hard functions = 0(a s ). We will not consider 
them in the following, since they are bound to be tiny. 

4.3 Analysis of the Q\ — Qj 7 contributions 

Evaluating the diagrams in Figure 4 for this pairing of operators, we obtain the direct photon 
contribution 

F$(Rr,ri = (-?) £ duS{u^) , (36) 

which involves a convolution of the leading shape function (5) with a jet function consisting 
of the cut of a hard-collinear loop. This is the expected extension of the first term in the 
expression for Ff 7 T (E^,fi) given in (7). 

Much more interesting is the single resolved photon contribution arising from this operator 
pair. As shown in the left graph of Figure 5, it is obtained by combining the 0(\ 3 ) SCET 
four-quark operator in (25), which contains two anti-hard-collinear quark fields in addition 
to a hard-collinear strange quark and a heavy quark, with the leading-power contribution 
(23) descending from Q 7l . According to the rules (21), the conversion of the two anti-hard- 
collinear fields into a photon and a soft gluon costs a factor of A 1 / 2 , giving a total suppression 
with respect to the leading term of A ~ Aq C d/^6- When the (anti-)hard-collinear fields 
are integrated out in the second matching step, we obtain the HQET diagram shown in the 
right graph of Figure 5. Here and below, horizontal (vertical) dotted lines represent Wilson 
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lines along the n (n) direction. The HQET matrix element corresponding to the graph on 
the right in Figure 5 contains a soft gluon field in addition to the two heavy quarks. In the 
language of the factorization formula (3) this corresponds to a single resolved contribution 
with a subleading shape function. For the contribution from the operator Q\ we obtain 

F$l(E„ /x) = \ (1 - 5 U ) [ A duj 5(u + p+) 

(37) 



ui + ie l \ 2E y ui J 



01 7 (w,Wi,Ji) , 



where the CKM-suppressed parameter 5 U has been introduced in (13), and we have defined 
the subleading shape function 



^e'^ (38) 

2,71 



(B\(hS n )(tn)fl(l + 7b)(^5„)(0)<^ {StgGfS n )(rn) (fiJ/i)(0)|S> 



2M B 

The penguin function F(x) has been defined in (9). For < x < 1/4 this function develops 
an imaginary part. Note that we adopt a power counting for the charm-quark mass such that 
m\ = C(m;,AQCD)- The argument of the penguin function in the convolution (37) then counts 
as 0(1). 

The structure of the soft Wilson lines in (38), which are directed either along n or n, follows 
when the decoupling transformation is applied to the (anti-)hard-collinear fields in SCET to 
remove their soft interactions from the effective Lagrangian and absorb them into eikonal 
factors. The Wilson lines reflect the space-time topology of the HQET diagrams shown on 
the right-hand side in Figure 5. Let us label the two weak vertices by coordinates (left) and 
x = tn + x + + x± (right), and the vertex of the soft gluon by y = rn + y_ + y±. The multipole 
expansion of the effective-theory fields implies that x +) ± and y~,± can be set to zero at this 
order. Gauge invariance then requires that the fields h(tn) and G s (rn) are joined by a Wilson 
line, and the rules of SCET determine that this line consists of two segments: a straight line 
[tn, 0] along the light-like direction n followed by a straight line [0, rn} along the light-like 
direction n. The fields G s {rn) and h(0) are joined by a straight Wilson line [rn, 0] along the 
light-like direction n. Using that [tn, 0] = S^tnjS^O) etc., we recover the structure of the 
Wilson lines in the non-local operator in (38). We note for completeness that soft functions 
closely related to our functions in (38) and g u in (58) were introduced, in a context not 
related to B — > X s ^ decay, in [66] . 

There is more to the space-time structure of the soft operator that is worth pointing out. 
Since hard-collinear fields in SCET carry large momentum components, the particles created 
by these fields always move forward in time. As a result, after convolution with the jet 
functions, the quantum fields in the definition of the subleading shape functions are ordered 
in the same way as they appear in Feynman graphs [32]. The operators considered in the 
present paper contain fields that propagate along the two light-like directions n and n, as 
indicated by the dotted lines in the right graph in Figure 5. If we assign coordinate to the 
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first of the weak vertices in the figure, the gluon is emitted at space-time point rn with r > 0. 
This is ensured by the ie prescription in the jet function in (37), since 



/ 



du\ e 



u>i + ie 



i 



2n9(r). 



(39) 



The gluon thus lives at a later time than the field h(0) (recall that n° = n° = +1), and indeed 
it appears to the left of that field. 

Another comment is in order concerning the structure of the result (37). From the dia- 
grams shown in Figure 5 we derive one half times the expression (37) without the real-part 
prescription in front of the integral and in expressions (4) and (13). The mirror diagrams 
not shown in the figure, in which the two weak vertices are interchanged, give an analogous 
contribution with the complex conjugate Wilson coefficients and CKM matrix elements, the 
complex conjugate penguin function F*(r), 2 the propagator factor l/(u>i — ie), and the soft 
function 



which is related to the original one by complex conjugation: g' 17 (uj, u>i, /i) = [g 17 (u, u>i, //)]*. To 
show this, one uses translational invariance to shift all position arguments by — tn and then 
changes the sign of the integration variable t. The sum of the diagrams in Figure 5 plus their 
mirror graphs thus gives a real result, and after averaging over CP-conjugate decay modes we 
obtain (37). 

The real-part symbols in (13) and (37) refer to different kinds of complex parameters. The 
various products of Wilson coefficients and CKM factors carry, in general, CP-violating weak 
phases. The convolution of the jet and soft functions, on the other hand, can carry CP-even, 
strong-rescattering phases, which in principle can result either from anti-hard-collinear loops 
(i.e., the jet functions Jj) or from the soft matrix elements themselves. (However, in Section 5 
we will argue that the soft functions are real.) When both types of phases are present, a non- 
zero direct CP asymmetry arises. The subleading power corrections investigated in this paper 
provide new mechanisms for generating such an asymmetry. A more detailed exploration of 
their phenomenological relevance is left for further study. 

The range of support of the soft functions in HQET can be derived by noting that the 
light-cone projections n ■ pi and n • Pi of all parton momenta in the B meson must be non- 
negative, and that the total momentum of all partons in the B meson is M B v. Since in HQET 
the momentum of a heavy quark is decomposed as pb = m^v + k, where k is the residual 
momentum, it follows that 
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(B\ (hS n ) (tn) {-irinp) (4 gGfS n ) (tn + rn) (S^S, 



n 



)(tn)fi(l + l5 )(Sth)(0)\B) 



2M B 




(41) 



2 This is because the fields in the charm-quark loop of the mirror graphs are anti-time-ordered. 
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where n ■ k > —m b and n ■ k > —m b . In the heavy-quark limit m b — > oo it follows that 
— oo < n ■ k < A and < n ■ pi < oo (for % ^ b), and similarly for n ■ k and n ■ pi. In 
the special case of the soft function g 17 in (38), the variable u corresponds to the residual- 
momentum component n • k of the initial-state heavy quark, while u\ can either correspond 
to the component n ■ p g of a gluon in the final-state B meson or to the component — n • p g of 
a gluon in the initial-state B meson. It thus follows that — oo < u < A and — oo < u>i < oo. 
This implies that the penguin-loop function F(x) is sampled over both positive and negative 
values of its argument. 

In principle, each of the six QCD penguin four-quark operators in the effective weak Hamil- 
tonian can give rise to a similar contribution, either via loops of massless quarks (q = u, d, s) 
or via a charm-quark loop. (6-quark loops lead to further power suppression.) Of these options 
only Ql has both a large Wilson coefficient C\ ~ 1 and a large CKM factor V cb V* s = 0(X 2 ), so 
it will give rise to the dominant effects. However, for academic reasons we will also consider 
the case of the CKM-suppressed operator Q± . Using that F(0) = 0, it follows that the resolved 
photon contribution resulting from this operator is given by 

2 C f°° did 

F^ u (E„fi) = -5 u Re du5(u+p + ) —i-g^u^fi . (42) 

6 J-oo J-oo UJ l+ l£ 

Note that the soft function is the same as in (37). 

Let us now investigate the convergence properties of the convolution integrals in (37) and 
(42). In the UV region, for oo\ 3> Aq CD , the first integral approaches the form of the second one, 
since mass effects become negligible. It follows that the convolution over io\ converges as long 
as the soft function g 17 vanishes for oo\ — > ±oo. In general, the asymptotic behavior of the soft 
functions for large values of the Ui variables can be analyzed using short-distance methods 
[29]. This shows, for instance, that the leading shape function behaves as S(oj,fi) ~ 1/ui 
modulo logarithms for u — > — oo. For the present case, naive dimensional analysis suggests 
the behavior g 17 (u, u>i, /i) oc u>i for large oj\ but fixed u, in which case the convolution integral 
would diverge linearly. To obtain such a contribution, however, would require a non-zero 
matrix element of the soft operator between two on-shell b quarks. But this matrix element 
vanishes by Lorentz invariance. A non-zero contribution is only obtained if, in addition to the 
heavy quarks, one adds a soft external gluon. This costs two orders in power counting, so that 
the asymptotic fall-off is at least as strong as gn(uj, u> 1 , ji) oc 1/Vl for u>i — > ±oo. It follows 
that the convolution integrals (37) and (42) are UV convergent. 

The behavior of the soft functions in the IR region cannot be derived from a perturbative 
analysis. In the present case, however, it suffices to make the reasonable assumption that 
gi 7 (co,Ui, fi) is non-singular at u>i — 0. Using the expansion 

1 - Fix) = l — l — - . . . (43) 

v ; I2x 90a; 2 560x 3 v ; 

valid for large x, we find that for small U\ the convolution integral (37) arising from the 
charm-quark loop behaves as 

r dui l _ p f m ie \ 17 ( w ~ __^2_ f dw 1 g 17 (u,uj 1 ,fj 1 ) . (44) 

J V 2 -^7 w i / J § m c J 
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For the convolution integral (42) arising from the up-quark loop, we find instead 



— g 17 {uj, ui, /I) = P / gn{w, UJ!,/!)- in g 17 {uj, 0, //) , 

u>i + is J Ui 



(45) 



where the symbol P denotes the Cauchy principal value of the integral. We conclude that the 
convolution integrals indeed exist as long as the subleading shape function is non-singular at 
Ui = 0. Note that for the case of the up-quark loop it is important that the integral over u> 1 
runs over both positive and negative values. Previous authors have already pointed out that 
the up-quark loop contribution to the B — > X s 7 decay rate, while strongly CKM suppressed, 
is described by an uncalculable long-distance contribution [35, 37, 38]. Our relation (45) 
provides a rigorous field-theoretic definition of this contribution in terms of a well-defined, 
non-local soft matrix element. 

For phenomenological purposes it is useful to define a new function 



/i7, ? (w,/x) 



u>i + ie 



1 — F 



IE 



(m h + ui) ui 



(46) 



Our final expression for the Q\ — Q 7l contribution can then be written as 



'17 



\ /J —p+ q=c,u 



(47) 



Note that the argument of the penguin function entering fi7, q (—p+, /i) is rn? q /[{mb — p+) wi] = 
m 2 q / '(2E 1 uji), as it should be. 

Our result differs from the parton-model expression given in (7) by the now familiar integral 
over the leading shape function in the first term, and by the fact that a non-trivial function 
f n, c (—p+, /-0 replaces (— 5(p+). The latter expression would be obtained if one made the 
unjustified approximation of neglecting the dependence of the soft fields in (38) on t and r, 
i.e., of evaluating all fields and Wilson lines at the origin. This would replace 



(48) 



and from (46) we would then recover the parton-model expression given above. However, there 
is no reason why (48) should provide a decent model for the soft function. All we know is that 
the definition of the soft function gn in (38) implies the normalization condition 



dcu duj igi7 {u,u u n) = ,a „f/ s ^ = 2A 2 , 

J— oo J— oo 



2Mf 



(49) 



where we have used a general relation derived in [67] to evaluate the matrix element of the local 
quark-gluon operator in terms of the hadronic parameter A 2 . Moreover, the trace formalism 
of HQET [44, 67] implies that the soft function can be written as 



g 17 (u,u u n) = Tr 



fl(l + 75) «7a ^tt^ Ea± ( u > n > w > w i> A*) 



4E 2 (uj,uj 1 ,fj,) , (50) 



23 



where we have used that the most general decomposition of the quantity E a± is of the form 
E a± (v , n, u, oji, p) = i / ~f a± ('E Jl +^S 2 ) with scalar functions S, = Sj(u;, a>i, p). It follows from 
this argument that the factor (1 + 75) in (38) can be replaced by 1, since the part of the trace 
involving 75 vanishes. It is then easy to see that 

dwg 17 {u,u u ij) = / 017(0;, p) . (51) 

—00 J —00 

One can constrain the function g 17 (u,Ui, p) further by looking at its first moments with 
respect to u and ui. These can be related to linear combinations of HQET matrix elements 
with three covariant derivatives. Such matrix elements can be expressed in terms of two 
hadronic parameters, p\ and p 2 , via [68] 

{B\hT a& ^DWh\B) _ 1 ^ /, i + y ^ _ ^ ^ * ^ p i] i + y 



2M £ 
Thus, we find 



a<5/3 



3 2 



(52) 



f K , f°° , , v (B\hii-f La in- D[iD1,in- D]h\B) 

/ rfww / duxgiriu^ufj.) = — = -p 2 

J-00 </-oo 



00 

A 



, ' (5|^7 ±a [[i^,m-D],m-D]/i|5) 
dw / dwiWi^i7(w,a;i,Ai) = ^ = 0, 



(53) 



where the HQET parameter p 2 is related to the parameter p\ s introduced in [69] via p\ s = 3p 2 . 
The vanishing of the first moment with respect to oo\ of g\ 7 is not a coincidence. As we will 
see in Section 5, g 17 is in fact a real function. Relation (51) then implies that all the odd 
moments in u> 1 vanish. 

As a final comment, let us add that even in the limit where the charm quark is treated as 
a heavy quark, m c = 0(rrib), the penguin contribution to the photon spectrum must still be 
described by a subleading shape function. In this limit the argument of the penguin function 
in (37) is of order m fe /AQ CD . Expanding then the function [1 — F(x)] to first order in \jx leads 
to 3 

m b + co fdt _ iwt (B\{hS n )(tn)j(l +75)^ {StgGfh)(0)\B) 

^'^"i^y 27 e 2M~ B • (54) 

Integrating this expression over lu, and dropping higher power corrections, we obtain a con- 
tribution to the total decay rate proportional to 



m b X, 



j _ dp+f 17jC (-p + ,fj.) ->■ (55) 



This agrees with the result found in [35-38] (the correct sign was obtained in the last reference) . 



3 Only the first term in this expansion gives rise to a UV-convergent convolution integral. 
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Figure 6: Examples of SCET diagrams giving rise to resolved photon contributions 
suppressed by at least two powers of of l/m&. The left graph arises from the pairing 
Ql — Q\, while the right one contributes to the Q\ — Q$ g term. 



4.4 Analysis of the Ql — Ql and Ql — Q 8g contributions 

The power-counting rules described in Appendix C show that for these two cases there do not 
exist operators arising at order 1/mi, in the heavy-quark expansion that contain soft fields 
other than the two heavy quarks. In particular, the diagrams shown in Figure 6 contribute to 
the B — > X s 7 photon spectrum only at order l/m^. This is an important finding. Since for the 
case of two charm-quark loops the first diagram in the figure is proportional to the large Wilson 
coefficient |Ci| 2 ~ 1, it has sometimes been mentioned as a potentially dangerous source of 
power corrections. It follows from our analysis that this contribution scales as (Aqcd/^6) 2 
relative to the leading term. It is therefore expected to be a small correction (see below). 

It thus remains to calculate the leading power corrections to the direct photon term in the 
factorization formula, by analyzing graphs of the type shown in Figure 4. Note that here only 
the first diagram on the left in this figure can contribute. After a straightforward calculation 
(summing over q = c,u), we find once again contributions involving a convolution of the 
leading shape function with a jet function consisting of the cut of a hard-collinear loop. The 
results are the same in the two cases and given by 



2E 1 



C F a s (n) 2 
4tt 9 



dui S(oj, /i) . 



(56) 



p+ 



This is the obvious generalization of the parton- model results in (7). Note that the prefac- 
tor mj,/2E 1 in the result for F 18 follows from the SCET calculation, and we have therefore 
presented it here. In the endpoint region this factor equals 1 up to power corrections, and it 
could therefore be omitted, as we have done in (12). 

In order to substantiate the statement about the smallness of the 0(1/ ml) double resolved 
photon contribution represented by the first diagram in Figure 6, we have evaluated this graph 
explicitly. The result is 



F 1 ( ? ) (E 7 ,/i) = - — 



1 



At 



duj S(u + p + ) / du\ 



(57) 
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g 11 (u,uj 1 ,U2,fi) , 
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where the l/m b prefactor indicates the additional power suppression. We have defined the 
soft function 



g 1 i(u,ui,u 2 ,fi) 



dt 
2^ 



-iuit 



\y^ 2U 



(58) 
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(B\(hS n )(tn) {SlgGfS n )(tn + un)T(SiS n )(tn) {SjS n )(0) (Sj gG^S^rn) (S^h)(0)\B) 

2M R 



where T = fl(l — 75). This function satisfies gu(ou, u>i, u>2, /j) = \gn(u), 002, u>i, //)]*, which 
implies that is real. In order to obtain an estimate of the magnitude of this contribution, 
we expand the penguin functions to first order using (43). This yields 



1 (2E- 



648 I m b 



Ac 
At 



dui \ dw 2 g u (-p+,Ui,U2,fJ>) , (59) 



where the remaining double integral over the soft function scales like Aq CD . For any reasonable 
value of this quantity, the prefactor 1/648 and the additional l/m b suppression render this 
contribution negligible. For instance, if we model the double integral by A n S(—p + ,fi) with 
some hadronic scale An ~ Aqcd, the contribution of this term relative to the leading direct 
photon contribution in (4) is approximately given by 



1 

648 



CM 

C 7y (n) 



A11 



-2 • 10" 



An V 
0.5 GeV/ 



(60) 



4.5 Analysis of the Q 8g — Q 8g contribution 

As we shall see, this contribution is more subtle than the remaining ones, so we will present 
its calculations in more detail. As mentioned in Section 4.1, the relevant matching relations 
for the dipole operator Qg g contains three SCET operators: the 0(X 3 ) four-particle operator 
in (26) containing the photon along with a hard-collinear strange quark and a hard-collinear 
gluon, and the two leading-order operators in (24) containing either a hard-collinear strange 
quark and an anti-hard-collinear gluon, or an anti-hard-collinear strange quark and a hard- 
collinear gluon (see the second row in Figure 3). After the conversion of the anti-hard-collinear 
partons into a photon plus soft particles only the second operator contributes at 0(X 3 ), while 
the third one receives a stronger power suppression. Note, in particular, that matching Q 8g 
onto an operator containing any soft fields does not give rise to a contribution at order l/m b 
in the heavy-quark expansion. It follows that the Q 8g — Q 8g term receives two contributions: a 
direct photon contribution from a pair of two SCET operators of the form shown in (26), and 
a double resolved photon contribution from a pair of two SCET operators of the form shown 
in the second line in (24), followed by the 0(\ 1 / 2 ) conversions of the anti-hard-collinear quark 
fields into photons. The resulting SCET diagrams are shown in the left panels in Figure 7. 
When the (anti-)hard-collinear fields are integrated out in the second matching step, we obtain 
the HQET diagrams shown in the right panels. 

The direct photon contribution involves a convolution of the leading shape function in (5) 
with a subleading jet function consisting of the cut of a hard-collinear loop. In the present case 
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Figure 7: Diagrams arising from the matching of the Qsg — Qsg contribution onto SCET 
(left) and HQET (right). Dashed lines denote non- localities obtained after (anti-)hard- 
collinear fields have been integrated out. 



the jet function is divergent and needs to be regularized. Using dimensional regularization 
and subtracting its 1/e pole in the MS scheme, we obtain 

P (a) /P \ C F a s {n) ( m b \ 2 f A (2 m b {u + p + ) 14 \ 

F ^ E -r^) = —^{2Ej J. p+ dUJ W n ^ + - Q -- Q c RS )S(u,»). (61) 

The scheme-dependent constant crs vanishes in the MS scheme, c^jg = 0. If, one the other 
hand, we adopt the dimensional reduction scheme, in which the Dirac algebra is performed in 
4 rather than d = 4 — 2e dimensions, then cp^ = 1. We will see later that the final answer for 
F 88 is scheme independent. 

The double resolved photon contribution gives rise to a more complicated structure, as the 
resulting soft matrix element contains four quark fields located at different space-time points. 
We find 

F « A = 9 ) J J" 5{U + P+) j ^ WTVs J „ W^Te 9 ™ {U > Wl ' " 2 ' ' 

(62) 

where we have defined the subleading shape function 

„cut/, ... , , ,,\ _ f dr -iu ir f du ■ f dt iujt 
g 88 (u^ui^) - J —e J —e J —e ^ 

jB\{hS n )(tn)T A {slSn) (tn) T n (sjs) (tn + un) \ X s ) (X s \ (sS n ) (rn)T n (sjs n ) (0) T A (sjh) (0)\B) 
X 2M B 
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As before, the Wilson lines render the soft matrix element gauge invariant. The sum over 
soft intermediate states X s with strangeness S — — 1 in the first equation arises since in this 
particular case the hard-collinear jet does not contain the strange quark. Note that only 
color-octet partonic states contribute to the sum, not physical hadronic ones. Performing the 
complete sum over states gives rise to the second equation, in which the two strange-quark 
fields are not time ordered but appear in the order shown in the formula. In the definition 
above 

r» = f(l + T5), r. (64) 

are projectors onto two-component light-quark spinors. In deriving the result (62) we have 
simplified the Dirac structure using the identity 



H^iL® 7±m ^ l±a = (d - 2f LL ® LL , (65) 

where d is the number of space-time dimensions. 

According to the discussion of the previous section, it follows that ggg* has support for 
— oo < uj < A and — oo < u>i j2 < oo. Note the difference in the sign of the ie terms in 
the two anti-hard-collinear propagators in (62), which is due to the fact that the anti-hard- 
collinear fields connected to the right weak vertex in the diagrams are anti-timed-ordered. 
Consequently, after convolution with the jet functions the position variables r and u in (63) 
are restricted to positive values, such that the fields in the soft (X s \ . . . \ B) matrix elements 
are time ordered, while those in the (B\ . . . \ X S ) matrix elements are anti-time-ordered, as it 
should be. Finally, we observe that the second equality in (63) implies the relation 

- * 

^f(a;,wi,W2,/i)J = g^(u},u}2, wi,fi) , (66) 

and since the convolution in (62) is symmetric in uji and u 2 up to complex conjugation, it 
follows that the final result is real. Contrary to (49) there is no useful normalization condition 
for the soft function g^. 

For phenomenological purposes, it will be convenient to define a new, real function 

x / s 2 f°° du! f°° du 2 cut/ 

/ 88 (w,/i) = -/ — / -g 88 {u,uj 1: uj2,h) , (67) 

9 J^oo + is J.oo oj 2 - %e 

in terms of which the second contribution to the photon spectrum is simply 

Fi\E„fi) = 47ra s (/i) (j^j /«,(-?+, . (68) 
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+ u) / — -V / — ^(w,^,^) (69) 



Note that the poles at u\ — and w 2 = are regularized by the prescriptions, in analogy 
with (45). 

Let us briefly comment on the structure of our result in light of the general factorization 
formula (3). The present case is our only example of a double resolved photon contribution 
at order l/m b . The jet function for the hard-collinear gluon is given by the cut of the gluon 
propagator, which up to trivial prefactors yields J{p 2 ) = 5{p 2 ) with p 2 = m b (u + p + ), in 
analogy with the tree- level expression for the quark jet function in (4). The jet functions for 
the two anti-hard-collinear quark propagators are, up to a trivial numerator factor, given by 
J{p 2 ) = l/(p 2 + ie), where p 2 = 2E 1 ui 1}2 in the present case. Hence, the triple convolution 
can be recast in the form (omitting scale dependences for brevity) 

= H J m b du J(m b (p + + to)) J 2E 1 d<jj x J(2E 1 Ui) J 2E 1 du 2 [J{2E 1 u) 2 )] * g%t(w, w 2 ) , 

in agreement with the factorization formula (3). The hard matching coefficient H = 1 at tree 
level. 

Our presentation above has hidden an important subtlety. The scale dependence should 
cancel (up to terms of order a 2 ) in the sum of the two contributions (61) and (62), but in 
order for this to happen the convolution (67) must contain a /z-dependent term at zeroth 
order in the strong coupling. This fact is incompatible with a multiplicative renormalization 
of SCET operators in the usual (convolution) sense. The resolution of this puzzle is that 
the convolution integrals over the soft function themselves are not convergent. In order to 
demonstrate this, we calculate the asymptotic behavior of the soft function for large values 
u>i t2 ^> Aqcd, corresponding to highly energetic light quarks. This behavior can be extracted 
using short-distance methods [29]. At leading order in perturbation theory, we simply need 
to replace the light-quark fields in the definition (63) by a cut propagator and perform some 
phase-space integrations. Working in d = 4 — 2e dimensions, we obtain 

cut/ \ d(u}i)u}} > /" A , , „, , \ / i \—e 

g 88 {u,u u u 2 ,n) = 7-: , 2 _ e ~-r r d(wi - u 2 ) / du S{u , n) (u - u) +.... 

wi,2»a QC d (47r) 2 e r(l-e) J w 

(70) 

Corrections to this result are suppressed by powers of a s or Aqcd/^1,2- The limit e — > 
is smooth and gives rise to a dependence ggg* oc ui6(u)i — ui 2 ). It is then obvious that the 
double convolution integral in (67) is logarithmically divergent in the UV region. 4 When the 
convolution is understood in the usual sense as an integral over renormalized functions, then 
this divergence is not regularized. 

On the other hand, our explicit expression (70) shows that the convolution integral would 
be regularized by the dimensional regulator if the limit e — > was taken after the convolutions 
have been evaluated. In that case we obtain a 1/e pole from the UV-divergent convolution 
integral, which needs to be subtracted in the MS scheme. We thus proceed as follows: we 



4 Note that according to (70) there are no UV divergences from the region of large negative values of Wi,2. 
In this region the convolution integrals are cut off by non-perturbative dynamics. 
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introduce a hard cutoff Auv an d split up the convolution integral in a low-momentum region 
defined by ooi,u 2 < Auv an d a high- momentum region defined by the complement. In the 
high-momentum region we can replace the soft function by the perturbative expression (70) up 
to higher-order terms in a s and power-suppressed contributions. We then evaluate the high- 
momentum contribution to the double convolution integral before taking the limit e — > 0. In 
doing so, we must remember to reinstate a factor (1 — e) 2 from the Dirac algebra, see (65), 
and a factor fi 2e from the conversion of the bare coupling constant g 2 into the renormalized 
coupling 4ira s (ii). In this way, we obtain 



/ 88 (w,/i) = - 



/-. \2 2e f°° f°° du 2 cut.bare/ 
v!-e) H / — / -^88 (U,U} 1 ,W 2 ) 

J-oo + *£■ J- 



/_ oc Ui + IE U 2 - IS 

'•Attv j, . /-Ai- 



ms subtracted 



2 f Avv dun r uv du 2 rut/ , , N 

- ^ jf V fj) (in Auv( ^-" } + 2 - 2c RS ) , 

where crs is the same scheme- dependent constant as in (61). This expression is independent of 
the auxiliary scale Auv, which for consistency should be taken to be several times Aq C d, so that 
perturbation theory can be trusted. In the above result the dependence on the factorization 
scale of dimensional regularization is explicit, and it is now evident that the sum of the two 
contributions (61) and (62) is both scale and scheme independent. Indeed, we find 



(72) 

8 ,( m A 2 dU ± dUJ 2 cut/ \ 

The large logarithm ln(mf,/Auv) i n the first term results from the ratio of the hard-collinear 
scale m b (u + p + ) in (61) and the soft (yet perturbative) scale Auv(^ + P+) contained in 
the function /s8( — i n (68). Resumming these large logarithms would require solving 

evolution equations in the effective theory. In the case of UV-divergent convolution integrals, 
the derivation of such equations is an open problem. 

Note that our result (72) is insensitive to the mass of the strange quark, as it should 
be. This is in contrast with the parton-model expression derived in [41] and shown in (7). 
The IR regulator m s introduced in the parton-model calculation is replaced in real QCD by a 
subleading shape function, i.e., by a hadronic matrix element of a non-local operator. We have 
checked that one recovers the parton-model expression for Fg$ if one calculates the soft matrix 
element in perturbation theory, i.e., if one assumes the validity of (70) also at small values of 
u>i >2 and introduces m s as an IR regulator, which replaces wf e (u/ — uj)~ £ — > [ujiiuj' — ui) —m 2 ]~ e 
in this formula. Of course, such a treatment cannot be justified due to the non-perturbative 
nature of QCD at low energies. 

It was argued in [34] that the IR-sensitive terms in the Q$ 9 — Qs g contribution to the 
B — > X s 7 photon spectrum can be absorbed into photon fragmentation functions of a strange 
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Figure 8: Diagrams arising from the matching of the Q^ — Q^g contribution onto SCET 
(left) and HQET (right). Dashed lines denote non-localities obtained after (anti-)hard- 
collinear fields have been integrated out. 



quark or gluon. While no formal proof of this assertion was given in that paper, it is likely to 
be true in the kinematic region away from the endpoint, where the splitting processes s — > 7 + s 
and g — > 7 + g can be treated using the collinear approximation. In the language of SCET 
this means that the partons after the splitting are still anti-hard-collinear fields, and hence the 
photon energy cannot be near the endpoint. In the endpoint region, on the other hand, these 
partons are soft, and they do not factorize from the remaining soft matrix element. Hence, 
in this region the non-perturbative physics is encoded in a complicated subleading four-quark 
shape function rather than a simpler fragmentation function. 

4.6 Analysis of the Q'j 1 — Q$ g contribution 

Evaluating the diagrams in Figure 4 for this operator pair, we find the direct photon contri- 
bution 

F 7 <«>(£ 7 , „) = °££M ^ j fju „) , (73) 

which generalizes the parton-model result in (7). The case of the Q77 — Qsg interference term 
is special in that, even though the parton-model expression for F^ rt {E 1) jj) does not indicate 
any problematic feature that would call for non-trivial soft contributions, there actually do 
exist some 0(l/m b ) effects that are described by subleading shape functions. Moreover, these 
effects remain non-local even for the total decay rate [39]. 

In order to study the resolved photon contributions, we must combine either one of the 
two SCET operators in (24) arising from the matching relation for Q Sg with the leading-order 
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operator in (23) descending from Q 7l (see also Figure 3). In both cases, the conversion of 
the anti-hard-collinear fields gives rise to one or more soft quark fields. The relevant SCET 
diagrams are depicted in the left panels in Figure 8, while the corresponding soft graphs 
resulting after the second matching step are shown in the right panels. In the first case, the 
second soft quark is generated by an insertion of a subleading term in the SCET Lagrangian. 
Evaluating the first contribution in detail, we find 



F is (#7> A*) = y ™ s (fi) R e J du5(u + p + ) J 
x [<7re(w, ui,w 2 ,n)- fif (w, wi, u 2 , //)] . 



s r du x r du 2 

7 J— oo J— oo wi + J-oo w 2 - is (74) 



The soft function g 78 arises when the hard-collinear strange-quark line in the left diagram 
in the first row of Figure 8 is cut, while the function originates from the cut through 
the hard-collinear gluon line. In this latter case the soft strange-quark line must also be cut. 
Specifically, we define the corresponding subleading shape functions as 

n (,, n n\- t ^ P~^ r f dU f dt p -iuA 

978 ( W)Wl)W2)/ ,)-j-e J -e J -e 

(75) 
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J 2tt J 2n J 2n y 1 

(B\ {hS n ) (tn) T A T n (Sis) ((t + u)n) [sS n ) (rn) T n (slS n ) (0) T A (S* n h) (Op) 
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where T n was introduced in (64), and T n is defined in the same way as T n but with n and 
n interchanged. One half of the contribution shown in (74), but without the real part pre- 
scription, arises from the original diagrams, while the mirror diagrams not shown in the figure 
give the complex conjugate of the above expressions. The two results combined give a real 
result, as indicated in (74). Note that there is no need to insert a time-ordering symbol in 
front of the light-quark fields in the definition of g 7 $ in (75), since after convolution with the 
jet functions the integration variables r and u are restricted to take positive values, and hence 
the light-quark fields have a space-like separation. In the second equation in (76), on the other 
hand, the fields are not time ordered because the non-local four-fermion operator arises upon 
performing a sum over intermediate states, as shown in the first equation. 
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Consider next the contribution shown in the second row of Figure 8. In this case the soft 
light-quark pair can carry any flavor. We obtain 



F$£(E^,ii) = A-Ka s (n) ^-Re / du5(uj + p + ) / du l I du 2 — — (77) 

^7 J-oc J-oo J-oo Ui -U 2 + IE 



X 



( — \— + — \ gf 8 \uj,uj 1 ,U2,fi) - ( — l — l — J ^(w.Ul.Ws,/!) 

\uji + ie co 2 —ieJ \tui + ie u 2 — iej 



where we have defined the subleading shape functions 

(1)/ \ f dr -iwir f du Au, 2 u f dt -iu,t 



X 



(B\(hS n ) (tn) (StS fl ) (0) T A fi(l + 75) (4^) (0) T £ , e q {qS n ) (rn) fl T A (4?) M |S) 
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^(c, Wl , c 2 , /x) = / ^ | ^ | ^ e-* (78) 



x 



(S| (hS n ) (tn) (& n S n ) (0) T A fl(l + 75) (Slh) (O)T^, e ? (qS n ) {rn) #75 T^ 4 (fijg) (un)\B) 
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where the sum extends over light quark flavors (q = u, d, s), and e q denote the quark electric 
charges in units of e. One half of the contribution shown in (77), but without the real part 
prescription, arises from the original diagrams, while the mirror diagrams not shown in the 
figure give the complex conjugate of the above expressions. 

In these definitions the light-quark fields are time-ordered, as indicated by the T symbols. 
That this is the appropriate ordering can be seen as follows. After convolution with the jet 
functions, for the terms containing the propagator l/(o>i + is) in the second line of (77) the 
integration variables r and u are restricted to the range r > u > 0. These terms correspond 
to the Feynman graph shown on the left in the second row of Figure 8, in which q(rn) should 
appear to the left of q(un). For the terms containing the propagator l/(ui 2 — is) the integration 
variables are restricted to the range u > r > 0. These terms correspond to the analogous 
Feynman graph with the opposite direction of the fermion arrow on the light-quark line, for 
which q(rn) should appear to the right of q(un). Hence, the proper ordering is indeed the 
ordering according to (light-cone) time. On the other hand, arguments along the lines discussed 
in [70] suggest that the time-ordering prescription is, in fact, not required for forward matrix 
elements and fields at light-like separation. We assume in what follows that the T symbol can 
be dropped in (78). 

Very little is known about the complicated four-quark shape-functions defined in (75), (76), 
and (78). Following the general arguments presented in Section 4.2, we conclude that the soft 
functions g?g and have support for —00 < oj < A and —00 < uji j2 < 00. However, in the 
case of we must require that oui i2 > 0. Note also the symmetry property 



J — c 



dou 



g^ 5 \uj,uj u uj 2 ,i2) = / rfw5fj 8 ' 5) (w,w 2 ,wi,/i) , (79) 

J — 00 
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which follows from the definitions of the soft functions in (78). 

The fact that in the case of g^ the operators involve light quarks of all flavors offers a 
strategy for modeling their matrix elements between 5-meson states. Unlike the case of the 
four-quark operators encountered for g^, g 78 , and t^f, here it is possible to (very roughly) 
estimate the matrix element by inserting the vacuum intermediate state between the two light- 
quark fields. The "vacuum-insertion approximation" (VIA) is used extensively in the study 
of local four-quark operator matrix elements, and we see no reason why it should work less 
accurately for non-local operators. Following [39] , we then obtain 



J d U ^8 5) (W,W1,W 2 ,/X) ^ a = -e sp ec (l - J^J 0+(-Wl,/i)0+(-W 2 ,/i) , 



(80) 



where e spec denotes the charge of the spectator quark inside the B meson, i.e., e spec = 2/3 
for B ± , and e spec = —1/3 for B° and B°. The quantity F(fi) is the HQET matrix element 
corresponding to the asymptotic value of the product }'b\/M b in the heavy-quark limit [71]. 
Finally, 0+(cu,/i) is the leading light-cone distribution amplitude of the B meson [72]. It is 
a real function with support for u > 0, which vanishes at u = and asymptotically falls off 
like 1/oj modulo logarithms [73]. Useful forms for this function have been derived based on 
QCD sum rules [72, 74, 75], the relativistic quark model [76], and model-independent moment 
relations obtained using the operator-product expansion [73, 77]. The support of 0+(cu,/i) 
implies that only negative values of oo\ and oo 2 give rise to non-zero contributions in (80), 
which is in accordance with the fact that q (q) describes an anti-quark in the initial (final) 
state. 

To conclude this analysis, we define the phenomenological functions 

(I) 4 [°° dui [°° dco 2 , , 

3 J -00 ^1 + is y_ 00 lu 2 - ie L 

POO pOO -j 

/re ) (w,^)= / duj i j du 2 — (81) 
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In the VIA, we obtain 



j — c 
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j — c 



VIA 




duRe/SW) _ = -e spcc — ^ 1 - —A TT7 - , (82) 



where \b = J^°duj 0+(w, fJ>)/w denotes the first inverse moment of the B- meson light-cone 
distribution amplitude [49]. In terms of the functions ffy U \ the direct and resolved photon 
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contributions to the B — y X s ^ photon spectrum can be summarized as 



CVa: s (/i) mf, 10 
4tt 2EZ Y 



•A 



-^78(^7, A*) 



du S(u, y) 



-p+ 



(83) 



+ Ana s (y) — - Re f^(-p+, y) + ft 8 >(-p+, y) 



5 Constraints from PT invariance 

In the expressions presented in the previous section there are four potential sources of com- 
plex phases: weak (CP-violating) phases from the CKM matrix elements and the Wilson 
coefficients, and strong (CP-conserving) phases from the new jet functions Jj and the various 
subleading shape functions. The CKM phases are non-zero only for the Q\ — Q 7l contribu- 
tions (with q = c,u), where they are suppressed by two powers of the Cabbibo angle. The 
Wilson coefficients are real in the Standard Model, although they can be complex in many of 
its extensions (see e.g. [78]). 

A unique property of the resolved photon contribution is that the new jet functions Jj 
are given in terms of full propagators (dressed by Wilson lines) and not by cut propagators. 
As a result, these functions are in general complex and give rise to strong phases. Since the 



relevant scale of the jet functions is y2£^ Aqcd 5 which is perturbative in the endpoint region, 
these strong phases are calculable in perturbation theory. The other potential source of strong 
phases are the soft functions, whose phases are of a non-perturbative nature. In studying 
the effects of the resolved photon contributions on the rate and CP asymmetry in B — y X s ^y 
decay, it is important to obtain some handle on these non-perturbative phases. In order to do 
so, we employ the invariance of strong-interaction matrix elements under parity (P) and time 
reversal (T). 

Under the combined transformation PT, a spinor field ip(x) transforms as PTip(x) PT = 
ApTip(-x), where in the Weyl representation of the Dirac matrices Apt = — 7°7 1 7 3 up 
to an irrelevant phase factor. The soft Wilson line S n (x) in (5) transforms into S n (— x). 5 
Finite-length Wilson lines, as they appear in the definitions of the soft functions, transform 
as PT[tn,0]PT = PT S n (tn) 5^(0) PT = S n (-tn) S't(O) = [-tn,0]. Finally, the external 
5-meson states transform as PT\B(v)) = —\B(v)). Also, because the time-reversal trans- 
formation is anti-linear, matrix elements get complex conjugated under application of PT. 
Consider now the definition of the soft function g 17 in (38). Using the fact that the position- 
space strong- interaction matrix element is PT invariant, we find that 



5 Strictly speaking the lower limit of integration is also changed from — oo to +oo, but this provides an 
equally valid definition of the same object. 
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where we have used that A PT ifLi^ A PT = fii^ and A PT 07 5 ry^ A PT = —^75^7^. However, 
we have already argued after (50) that the term containing 75 vanishes. Hence, PT transforms 
the position-space matrix element into the complex conjugate of the same matrix element with 
all position arguments Xi replaced by —Xi. By taking the complex conjugate of relation (84) 
and reversing the sign of the integration variables r and t, it then follows that gn{oj, 001, //) is 
a real function. 

An analogous argument can be presented for the soft functions gj$ in (75) and gfy 5 ^ in 
(78), where it is important however that we can avoid the time-ordering prescription for the 
soft light-quark fields. The HQET trace formalism can be used to show that in the definitions 
of these matrix elements only even numbers of 75 matrices can give rise to non-vanishing 
contributions, and then the Dirac structures are in all cases even under PT. More specifically, 
in analogy to (50) we can write 



g 78 (u,u 1 ,u 2 ,fi) = Tr 



— ^— T A n) Y B — ^— ^ 2 {v, n) 



+ Tr 



7s T A ^ 3 (v,n) 



Tr 



^4{v,n) 75 



(85) 



where T A = T n and Tb = T n , and for brevity we have suppressed the dependence of the 

coefficient functions Sj on 00, 001, 002, and u. A similar expression, but with different matrices 

(1 5 s ) 

T A and Tb, holds for g 7 % after a Fierz transformation. The most general Lorentz-invariant 
decompositions of the functions involve products of up to four 0, 0, and 7° matrices, where 
all transverse indices must be contracted. No 75 matrices appear in this decomposition. Note 
also that the relation n + n = 2v allows us to eliminate in favor of and 0. With only two 
independent external vectors v and n, however, it is impossible to saturate the four indices of 
an e Q( g 7 5 symbol, and hence only even numbers of 75 matrices in the product structure T a ®Tb 
can give rise to non-zero traces. For the case of g 78 considered above, it follows that we can 
replace 16 T n ®T n — > 00 (g) 00 — 0075 <8> 007s- Both of these product structures are even 
under PT. In the case of g$ and g 7 f , we find similarly that 0(1 + 75) ® — > ® and 
0(1 + 75) ® 075 — >■ 075 ® 075. Once again, these product structures are even under PT. It 
follows that the functions g 7 $, g$ , and gr^ are all real. 

Let us finally consider the functions g^ in (63) and g^ in (76), which are defined in 
terms of sums over intermediate states \X S ), which without loss of generality can be chosen 
to be eigenstates of PT with eigenvalues ±1. After summing over the polarizations of the 
intermediate states and integrating over their momenta, we find that each term in the sum 
over states can be written as a product of two traces, in analogy to the second term in (85). 
The same arguments as above then show that g^ and g 7 ^ 1 are real. 

In conclusion, we find that all of the subleading shape functions are real. The strong 
phases mentioned in the introduction to this section thus arise only from the new jet func- 
tions Jj. Given that the soft functions are real, it now follows from (51), (66), and (79) that 
/ dou g 17(00,00 1, /j,) is an even function ofo>i,and that g^(oj, 001, 002, //) and J doo g 7 g 5 \oo, 001, 002, n) 
are symmetric under the exchange of oo\ and 002- 
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6 Partially integrated decay rate 



In the various parts of Section 4, we have derived explicit expressions for the direct and resolved 
photon contributions to the coefficient functions F i j(E 11 /x) entering the master formula for the 
B — > X s 7 photon spectrum in (4). The results are given in (33), (47), (56), (72), and (83). 
For phenomenological purposes, it is most interesting to study the partial B — > X s ^ decay 
rate 

F(E ) = dE, — , (86) 

Je 



dE^ 



obtained by integrating the photon spectrum over a region E < E 1 < Mb/ 2. Provided that 
A = rrib — 2E is much larger than Aqcd, the direct photon contributions to this integrated rate 
can be calculated in terms of local operator matrix elements [26] using a combined expansion 
in powers of A/m b and Aqcd/ A- In the limit E — )• one obtains the total decay rate, and 
A = m b ; however, the rates measured experimentally are obtained with values of E larger 
than 1.7 GeV, so that A < 1.25 GeV. 

An important feature of the resolved photon contributions studied in this work is that 
they do not reduce to local operator matrix element in the limit A ^> Aqcd- Rather, the 
corresponding contributions to the integrated decay rate must still be described in terms 
of matrix elements of non-local operators. This implies that the corresponding theoretical 
uncertainties do not reduce significantly as the cutoff E is taken out of the endpoint region. 
We will now illustrate this by deriving expressions for the first-order power corrections to the 
integrated decay rate 



r(£ ) 



G 2 F a\V tb V t * s 
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(87) 



valid for A ^> Aqcd- Here rrib denotes the pole mass of the b quark. The dots represent terms 
of order 1/m^ and higher, which we ignore. The integrated coefficient functions are obtained 

as 



i^(A,/i) 



-A 



dp + Fij(E^ii) , 



where 



2Ej. As will be explained below, with the exception of ggg* the non 



perturbative soft functions have support for values u = (9(Aqcd)- In the limit A ^> Aqcd, 
the to integrals in the definitions of the subleading shape function can then be performed over 
the entire range from — oo to A, and this leads to simplifications. However, the integrals over 
the remaining u>i variables cannot be simplified. 



6 We ignore radiative tails of these functions, which can exhibit power behavior and extend to larger ui 
values. These effects only contribute at higher orders in a s . 
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For the direct photon contributions, we need the integrals 



dui S(ui, fx) « A , 



p+ 



dpi 



-A 



A ^ln m ^+ P+) g^^Afln^ 

-p + 



(89) 



where the 
(A Q cd/A) 
to (87) at 
expansion. 

Let us 
integrated 
tributions 
From (31) 



approximate expressions on the right are valid up to corrections of order a s (A) and 
2 , both of which are known [26]. It follows that the direct photon terms contribute 
order A/mb in power counting and can be computed using a local operator-product 
. In the formal limit Eq — » these terms are promoted to 0(1) contributions, 
now discuss what happens to the subleading shape-function contributions to the 
decay rate. For the operator pair Q 7l — Q 7l , the subleading shape-function con- 
also reduce to matrix elements of local operators, as discussed in detail in [30-33]. 
-(35) it follows that 
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Note that, at order l/m 6 , the non-zero strange-quark mass effect discussed in Section 4.2 
integrates to zero in the partially integrated decay rate (at tree level in a s ), as long as A > 
Aqcd- Similarly, for the operator pairs Q\ — Q\ and Ql — Qs g , only direct photon contributions 
contribute at order l/m b , and we obtain 



(91) 



The remaining contributions, all of which contain resolved photon terms, are more inter- 
esting. For the operator pairs Q\ — Q 7l , we obtain 



where 
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The integral over p + in (88) eliminates the S(u + p+) distribution in (37), and integrating 
the soft function g 17 (ui,uji, /i) in (38) over oj then eliminates the ^-integral and sets t — 0, 
so that part of the non-localities of the operator are eliminated. However, the gluon field is 
still smeared out on the n light-cone. Note that there is no contribution from the up-quark 
penguin loop to the integrated rate. As noted in Section 5, the integral over u of gn(ou, u>i, /i) 
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is symmetric in oo\, so that the integral over F$ in (42) vanishes. 7 In the approximation 
where the penguin function is expanded to first order using (43), one would obtain 

A,(a„ )m £^(-!)a-(i-^, (m) 

which equals the integral over the partonic expression in (7), where we had neglected the small 
correction proportional to 5 U . 

For the case of the pair Q 7l — Q 8g , we obtain 

F 78 (A^) = ^M^A + 4 7 r^)Re f° f° -**L- fflfa,^) , (95) 

where in analogy with (93) we have introduced 

4sVi,^) = J Y^ r \ ( 96 ) 

(B | (fegfl) (0) T A # 75 (0) g g e q (qSn) (rw) # 7 s T A (4g) W | g) 
X 2M B 

Note that the contribution from g^ 8 vanishes, since the integral over u of this function is 
symmetric under the exchange of uj\ and oo 2 . Likewise, the contributions from the functions 
g 78 and to the integrated decay rate cancel each other. This follows from the fact that 
the two non-local operators in (75) and (76) coincide for t — 0. In the VIA we obtain 

F78(A ' via = " IT" T A " —2~ 6spcc v 1 - N* ) A|R • (97) 

The second term coincides with the result derived first in [39]. 

Finally, for the case of the operator pair Q 8g — Q 8g , we find from (72) 
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where 

r-A 



h^(A,u 1 ,u 2 ,fi)= f dugX{u,u u u 2 ,fi). (99) 

J-A 



7 It has been pointed out in the past that up-quark penguin loops might give rise to an C(Aqcd/to(,) 
uncertainty in the integrated rate for B — > Xdj decay [38], where unlike in B — > X s -y they are not CKM 
suppressed. Applying our analysis to B — > X<j7 shows that this contribution actually vanishes, removing that 
source of uncertainty in the integrated decay rate. Note that the same is not true for the CP asymmetry in 
B — > Xdf decay, where the corresponding contribution is proportional to hi?(0), which is non-zero in general. 
We further comment on B — > Xd"f decay in the conclusions. 
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Naively, we would expect that for A ^> Aqcd this function becomes independent of A and 
reduces to the expression 

^V,^) = J ^e-^| *t e *** (100) 

(B\ (hS n ) (0) T A (SjSn) (0) Tn {Sis) (un) (sS n ) (rn)T n {S\S n ) (0) T A (Sjh) (Op) 
X 2M B 

in which case the second term in (98) would be strictly positive. However, in the present 
case the limit t — > in (63) is singular, since then the separation between the two light-quark 
fields s and s becomes light-like. As a result, the integral over u in (99) diverges linearly 
as A is raised to infinity, and hence it must be evaluated at large but finite A. In other 
words, unlike for the other soft functions, the support of the function g^ 1 is not restricted to 
values oj = (9(Aq CD ) but extends to large negative values of u. This is in accordance with the 
asymptotic behavior derived in (70). 

The convolutions of the soft functions with anti-hard-collinear jet functions in the results 
given above cannot be expressed in terms of local operator matrix elements, but rather define 
unknown hadronic parameters of order Aqcd- These are the sources of genuine, first-order 
power corrections to the integrated decay rate, which are not reduced by lowering the cutoff 
Eq on the photon energy. 



7 Phenomenological implications 

The results of the previous sections can be used to quantify the effect of the resolved photon 
terms on the B — > X s ^y photon spectrum, as well as on the decay rate and CP asymmetry. In 
this paper we will restrict our attention to the decay rate integrated over a sufficiently wide 
energy range. The photon spectrum and CP asymmetry will be studied in a future publication. 

In order to estimate the irreducible theoretical uncertainty from these new non-local effects 
on the integrated decay rate, we define the function 

r E (A) = r^-rWIopE (1M) 

1 1-^0 J I OPE 

where E is the lower cutoff on the photon energy, and A = — 2E . This definition is such 
that the true decay rate T(E ) is obtained from the theoretical expression r(_E )|oPE obtained 
using a local operator product expansion by multiplying it with [1 + ^(A)]. Note that 
T(i?o)|oPE refers to the formula used in previous calculations of the B — > X s ^ rate, see e.g. [5]. 
The function J-"e(A) corresponds to the relative theoretical error made in these calculations 
due to the neglect of non-local power corrections from resolved photon contributions. 
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To the order we are working, we obtain 
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where we have assumed that the Wilson coefficients are real (like in the Standard Model) 
and neglected effects proportional to V ub . Note that the terms in the first line on the right- 
hand side vanish due to the relation (90). We can express the various other contributions in 
terms of suitably defined hadronic parameters of order Aqcd, using the expressions for the 
quantities Fij(A,fi) derived in the previous section under the assumption that A ^> Aqcd- 
Making explicit the dependence on the Wilson coefficients and factors of the strong coupling 
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In the case of A i7 and A 88 we have factored out the appropriate powers of the quark electric 
charges. Because of the sum over light-quark flavors in (96), the parameter A 7 g ec receives 
contributions proportional to any one of the light-quark charges. The resulting hard breaking 
of isospin symmetry implies that its value will be different for charged and neutral B mesons, 
even in the limit of exact isospin symmetry of the strong interaction. We will show in Sec- 
tion 7.2 that, in certain approximation schemes, A 7 | ec is proportional to the electric charge of 
the spectator quark in the B meson. 

Note that the parameters m 2 c /m b and A entering the arguments of A 17 and A 88 count as 
C(Aqcd)- The dependence on the strange-quark mass in (103) arises only because the function 
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J-e(A) is defined as the deviation from the partonic rate r part (_E ). The true decay rate T(E ) 
in (101) is independent of m s . Note also that the result for A 88 is formally independent 
of the UV cutoff Auv, and that it is the only hadronic parameter in (103) that depends 
on the quantity A. In the formal limit where the cut on the photon energy is removed, 
A m&, the linear growth (modulo logarithms) of the parameter A§8 with A implies that the 
corresponding contribution to J-^A) is promoted from a power-suppressed to a leading-order 
effect. Indeed, it is well known that in this limit there exists a leading-power, non-perturbative 
Q 8g — Q Sg contribution related to the photon fragmentation off a strange quark or gluon [34] . 
For practical applications this observation is irrelevant. We will argue in Section 7.3 that, for 
realistic values of E Q outside the endpoint region, the dependence of A 8 s on A is very weak, 
and therefore the function .Fe(A) is almost equal to a constant. 

Without further information about the soft functions, the parameters are expected to 
be of order Aqcd apart from the electric charges factored out in (104). This would lead to 
very large effects of up to 30% on the decay rate. Fortunately, it is possible to constrain the 
values of Ai 7 and A^ 60 by means of simple considerations, as we will now discuss. The input 
parameters used for the estimates in the following discussion are collected in Appendix B. 
The accuracy of our calculations is such that we are insensitive to the scale dependence of the 
subleading soft functions and the corresponding hadronic parameters. Even though we have 
indicated their \x dependence in the formulae given above, to properly control this dependence 
would require to extend our calculations to the next order in the expansion in powers of a s (/i). 



7.1 Analysis of the Q\ — Qi 1 contribution 

In order to obtain a reasonable estimate for the parameter A17, we first collect everything we 
know about the function hyj(uj\,ii) defined below (93). As proved in Section 6, this function 
must be real, and the symmetry relation (51) then implies that it is an even function of u\. 
It follows that all odd moments of h 17 vanish. Moreover, from (49) the normalization of hn 
is fixed to 2X 2 . About the higher even moments nothing definite is known, but we can expect 
them to be proportional to an appropriate power of Aqcd times a not too large numerical 
factor. Finally, as a soft function, hi 7 should not have any significant structures, such as 
peaks or zeros, outside the hadronic energy range. 

The first functions that come to mind are an exponential and a Gaussian, 

hiriuufi) = — , or h 17 (ui,fi)= , (105) 

o y/2-Ka 

for which all even moments are finite. As long as a < Aml/rrib ~ 1.1 GeV, which with the 
power counting adopted in this paper is formally of order Aqcd, then for all relevant ui\ values 
the argument of the penguin function F(x) entering the definition of A i7 in (104) is much 
larger than 1/4, which is the radius of convergence for the Taylor expansion given in (43). It 
is then a good approximation to expand the penguin function [1 — F(x)} to 0(l/x 3 ). The first 
term in this expansion corresponds to the non-perturbative correction identified in [35] , which 
was already included in the partonic result and subtracted in (104). It therefore does not 
contribute to J-£;(A). The next term gives rise to an odd moment of hyj and thus vanishes. 
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The third term in the expansion contributes the amount 



A e r nded = _|_^| A2( ^ ) (1Q6) 

to A 17 . Here (ouf) denotes the (normalized) variance of the function hn(ui,fi), which equals 
2cr 2 for the exponential form and a 2 for the Gaussian. For a typical hadronic scale o = 0.5 GeV 
this gives 

A ex P anded = _ 6 9MeV and _ 3 .4MeV, respectively. Here and below we have used 
the input parameters collected in Appendix B. The corresponding contributions to the decay 
rate are very small, below 0.5% in magnitude. 

It is interesting that, due to a weaker numerical suppression, certain 1/rrib corrections to 
A17 can give a contribution of comparable size. They arise from the fact that the first moment 
of the function gnioj, u>i, /i) with respect to u> does not vanish, see (53). In order to calculate 
the resulting power-suppressed term, we replace the first relation in (104) by 
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where the factor (V^Lf appears because of the prefactor in (4). Expanding now the 
penguin function to first order yields 

(108) 

where the dots represent higher-order terms in the expansion of the penguin function, which 
in particular give rise to the contribution (106). The expression shown above yields a l/m b - 
suppressed contribution to the parameter A17, which we denote by 8 An. It is proportional to 
the normalized first moment of the function gn with respect to u, which according to (49) 
and (53) is given by (u) = -p|s/(6A 2 ) ~ 0.24 GeV. We obtain 

6Al7 = Thrf ~ -(9.8±5.2)MeV, (109) 

which is formally a power correction proportional to Aq CD /m6 to the result in (106). Here 
Pls = corresponds to the spin-orbit term of the HQET Lagrangian introduced in (53). 

In practice, it turns out that (106) provides a reasonable approximation only as long 
as a < 0.3 GeV. Performing the convolution integral in (104) exactly, we find that for both 
model functions in (105) the resulting value of | A17 1 is maximized for certain values of a, which 
depend on the functional form of hn- Using the input parameters collected in Appendix B, we 
obtain (Ai7 P ) max = -4.6 MeV for a = 0.51 GeV with the exponential model, and (Af 7 auss ) max = 
— 8.1MeV for o = 0.77 GeV with the Gaussian model. Note that the maximum values are 
smaller in magnitude than those one would derive from (106) with these values of a. 
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The above estimates do not provide a conservative bound on the size of the hadronic 
parameter A 17 . A significantly larger effect can be obtained if the soft function g 17 (u,Ui, fi) 
exhibits a tail outside the region |o>i| <C Am 2 c j m^. In analogy with the leading-order shape 
function, we expect that the function #17 exhibits a radiative tail proportional to 1/oui for large 
u>i. But even at the non-perturbative level, it is conceivable that a significant contribution to 
the integral results from the region of larger oj\ values. Consider, as an example, the model 

2A 2 uj - A 2 _4 
M^H^^^e ^, (110) 

which for A and a of order Aqcd satisfies all requirements one would reasonably impose on the 
soft function. The solid curve in Figure 9 shows this function evaluated with a = 0.5 GeV and 
A = 0.425 GeV. It features regions of positive and negative values and hence is less constrained 
at larger u> 1 by the fact that the normalization is fixed to 2A 2 . Having values of either sign is not 
problematic, because there is no probabilistic interpretation of the subleading soft functions. 
The long-dashed line in the figure shows the weight function under the convolution integral in 
the definition of A17 in (104), including the charge factor e c . With the above parameter choices 
for the soft function, we obtain A 17 = —42 MeV. By using another set of values, a correction 
with the opposite sign and of the same magnitude can be obtained. For example, taking 
a = 0.5 GeV and A = 0.575 GeV we find A 17 = 27 MeV. If we include the 1/mt, corrections 
as shown in (107), using (m^ + u) — > (mi, + (00)) = (nib — p|,s/6A 2 ), we find —62 MeV and 
21 MeV, respectively. Of course, these are just illustrative values, and one could obtain even 
larger negative or positive values by reducing the separation between a and A, which however 
will also increase the value of the soft function at uoi = 0. Nevertheless, based on these 
considerations, it seems to us that 

-60 MeV < A 17 < 25 MeV (111) 

is a reasonably conservative range, which we will adopt for our analysis below. While this 
allows for a value significantly larger in magnitude than the naive estimate (106), it nev- 
ertheless strongly suggests that A i7 is considerably smaller in magnitude than Aqcd- Note 
that the effect of a value of A 17 near the extreme values indicated above would be of the 
same magnitude as the effect of the leading-order, non-perturbative correction [35] resulting 
from the term proportional to A 2 in the expression for Ffj Tt in (7), which corresponds to 
-m b X 2 /(9m 2 c ) « -48 MeV. 

7.2 Analysis of the Qi 1 — Q 8g contribution 

It is instructive to analyze this contribution using the language of flavor symmetry of the 
strong interaction. Due to the weighting by the quark electric charges, the relevant four-quark 
operator in (78) is a pure SU(3) octet, which can be decomposed into two parts corresponding 
to isospin / = 0, 1. The Wigner-Eckart theorem implies that 

AS" = I Affl, ± \ Afl ' = I ( A « " A '-) + A '- ' (112) 
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Figure 9: Model function hn(ui,fi) from (110) in units of GeV, with o = 0.5 GeV and 
A = 0.425 GeV (solid line); weight function under the convolution integral in the definition of 
A 17 in (104) in units of GeV -1 (long-dashed line); weight function including l/m b corrections, 
obtained by the substitution lo — > (u) in (107) (short-dashed line). See text for explanations. 



where the upper (lower) sign in the first equation refers to charged (neutral) B mesons, and 
as before e spe c denotes the electric charge of the spectator quark in units of e. In the limit 
of unbroken S77(3) flavor symmetry, it follows that A^ = A^ 1; since both parameters arise 
from the matrix element of the same SU(3) octet operator. Hence, in this limit we obtain 

A 78 \su(3) = es P cc ^=1 • (113) 

Interestingly, the VIA discussed in Section 4.6 also predicts that A s ^ ec is proportional to 
e sp ec [39] , and we can use this fact to obtain a model estimate of the relevant 577 (3) reduced 
matrix element. From (82), we read off 

A 5=iIvia = A ?2oU = " (l " ^ ) ^) e [-386 MeV, -35MeV] , (114) 

where in the last step we have used the parameter ranges discussed in Appendix B. 

According to (112), the isospin-averaged decay rate [T(B° — >■ X s ^) + T(B~ — >■ X s j)]/2 
depends only on A® , while the isospin difference [r(JB° — >■ X s ^) — T(B~ — > X s ^y)] is propor- 

(8) 

tional to A}^. While a priori these two non-perturbative parameters are unrelated, we have 
just shown that they coincide both in the 577(3) flavor-symmetry limit and in the VIA. It was 
first pointed out in [79] that, in the limit of exact S77(3) flavor symmetry, the isospin-averaged 
decay rate can be related to the isospin asymmetry, 

T(B° -» X sl ) - T(B- -» X sl ) 
°- T(B° -> X sl ) + T(B- -> X sl ) ' 1 1 
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without employing the VIA. This asymmetry has been measured by the BaBar Collaboration 
using two different experimental methods. For the "sum-over-exclusive-modes method" with 
£ 7 > 1.9 GeV, they find A _ = (-0.6 ± 5.8 ± 0.9 ± 2.4)% [23], where the errors are statistical, 
systematic, and due to the production ratio B°/B~ , respectively. For the "recoil method" with 
E 7 > 2.2 GeV, they obtain instead A _ = (—6 ± 15 ± 7)% [80], where the errors are statistical 
and systematic, respectively. The naive average of these two results is A _ = (—1.3 ± 5.9)%. 
To the order we are working, the parameter 1 is related to Aq_ via 

AflA =-^44-^-A _^(59±268)MeV, (116) 
/ " 1|exp C 8s (aO 27ra s (/i) 

where in the last step we have used the average experimental result with its large uncertainty 
given above. This value is consistent with the prediction (114) obtained in the VIA within 
errors, even though the central value has the opposite sign. 

Allowing for SU(3) flavor-symmetry breaking at the level of 30%, i.e. Ag = (1±0.3) Ag x , 
we finally obtain 

Al" = (e S pcc ± 0.05) Afl, « -4.5 GeV (e spcc ± 0.05) A _ , (117) 

which is meant as a range, not an error bar. This formula implies that, within the quoted 
uncertainty, the isospin asymmetry also determines the flavor-averaged value of A s ^ cc . For the 
corresponding contribution to the flavor- averaged value of the function ^(A), we obtain 



^ VS (A)L = -(1±0.3)^, (118) 



which adds S*?7(3)-breaking effects to the estimate derived in [79]. Note that this relation is 
independent of the values of the Wilson coefficients and other theoretical parameters. 

Due to the current large experimental uncertainties in the measurement of the isospin 
asymmetry, it is difficult to give a reliable estimate for A s ^ cc . Based on (114) and (116), we 
expect that the parameter Ag x is negative (assuming that the VIA is sufficiently reliable to 
predict the sign correctly), but since the experimental value allows for the entire range in (114) 
at the level of two standard deviations, we cannot restrict that range further at present. A 
future, more accurate measurement of A _ could improve the situation. 



7.3 Analysis of the Qs g — Qsg contribution 

Unfortunately, we have very little useful information about the soft function entering the 
definition of the hadronic parameter A 88 in (104). Its asymptotic behavior for large values of ui-y 
and Cc>2 can be derived from (70), and it ensures that A 88 is independent of the UV cutoff Auv- 
Note that the second term in the definition of A 88 , which contains the logarithm of Auv/ A, is 
bound to give a very small contribution to A 88 , because (CVe^A)/(87r 2 ) < 3MeV is very small 
for realistic values E > 1.6 GeV. We thus expect that the hadronic parameter A 88 receives 
its dominant contributions from values iO\^ = (9(Aqcd), for which no useful constraints on 
the soft function h^f exist. For the same reason, we expect that the linear growth of A 88 for 
large A is a numerically irrelevant effect. It then follows that the function ^^(A, u>i, u>2, //) is 
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approximately equal to the function h^(uji, u 2 , //) shown in (100), even though this relation 
is not strictly valid. As mentioned earlier in the paragraph following that equation, this form 
would imply that the contribution to A 88 resulting from the double integral in (104) were 
strictly positive. 

In summary, we expect that the hadronic parameter A 88 (A, fj) is, to a good approximation, 
independent of A and given by a positive, non-perturbative constant of order e 2 s Aq CD : 

A 88 (A,/i)^A(/i), A(/x)>0. (119) 

We have backed up this expectation by using different models for the soft function /igf, for 
example by writing it as a product of two functions f\{oJ\) 72(^2) and using various models 
such as exponentials or Gaussians. A particularly simple example is provided by functions 
hg$(ui, u>2, fi) that are symmetric in both ui variables and have support for oj-i = (9(Aqcd)- 
In this case the third relation in (104) implies A(fi) 27r 2 /i 88 t (0, 0, /i), and the value at the 
origin scales like /^(O, 0, fi) ~ Aqcd- For our numerical analysis, we will consider the rather 
generous range < A(/i) < 1 GeV. Even for the largest value, the suppression by the charge 
factor e\ = 1/9 in (119) renders the effect of this term on the decay rate to be very small. 



7.4 Summary of phenomenological estimates 

We are now in a position to study the implications of our analysis for the function J-'e(A) in 
(103). Using the parameter values collected in Appendix B, we obtain from (111) and (119) 
the contributions 

(120) 

F E \ 88 E [-0.3, +1.9]%. 

The value of J^eUs depends slightly on A and is obtained using A = 1.45 GeV, corresponding 
to a cut at E = 1.6 GeV. For the case of J-^ks; we consider the charge-averaged contribution 
and quote separately the theoretical estimate obtained using the VIA and the experimental 
estimate derived from the measurement of the isospin asymmetry. In the latter case we allow 
for 30% SU(3) violation, as indicated in (118), and take the 95% confidence level experimental 
range. This yields 

^|!J A G [-2.8,-0.3]%. 



(121) 

F E \ e ^ e [-4.4, +5.6] % (95% CL) . 

In order to obtain a conservative estimate of the combined theoretical analysis, we adopt 
a Baysean approach and add up the various contributions using the scanning method. In this 
way, we arrive at our final result 

-4.8% < F E (A) < +5.6% (VIA for A s 7 n , (122) 

where we have used the theoretical estimate for J^e^s- When the experimental estimate is 
used instead, the range is expanded to 

-6.4% < F E (A) < +11.5% (A s 7 P ec from A _) . (123) 
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We emphasize that the estimates in this sections should be considered as ranges, within which 
we expect the actual values of Te to lie, without making a statement about the most likely 
values within these ranges. 

If in the future a more precise value of the isospin asymmetry can be measured, this could 
be used to reduce the uncertainty range somewhat. If, for example, we assume that the true 
isospin asymmetry lies in the center of the interval predicted by the VIA, A _ = +4.6%, then 
in the absence of experimental uncertainties we would derive J-'eIts' £ [— 2.0, —1.1] %, where 
the remaining uncertainty stems from the unknown effects of SU (3) breaking. In this "ideal" 
case, the combined result would be 

-4.0% < T'b(A) < +4.8% (ideal case) . (124) 

We do not see a possibility to reduce this uncertainty in the foreseeable future, given that no 
theoretical tools exist to constrain the non-local matrix elements defining the soft functions 
entering the various resolved photon contributions studied in this paper. We therefore consider 
the range in (124) as the irreducible theoretical uncertainty affecting any theoretical prediction 
of the B — > X s 7 branching ratio. 

8 Conclusions 

The inclusive radiative decay B — > X s ^ is used extensively in constraining extensions of the 
Standard Model. For example, it provides very stringent constraints on extended Higgs sectors 
in type-II 2-Higgs-doublet models and supersymmetric models. The theoretical prediction for 
the corresponding branching ratio is at a stage of precision where the remaining perturbative 
uncertainties are estimated to be of order 3% [5]. The limiting theoretical uncertainty arises 
from non-perturbative effects outside the realm of the local operator product expansion [39]. 
It is therefore important to analyze these effects in a systematic fashion. In this paper, we have 
for the first time provided a complete analysis of non-local l/m& corrections to the B — > X s 7 
photon spectrum and decay rate, working at tree level in perturbation theory. Compared to 
inclusive semileptonic B decays, non-perturbative effects in radiative decays are much more 
complicated to analyze. First of all, one must consider the contributions of many different 
operators in the effective weak Hamiltonian, not just one operator. More importantly, however, 
new types of non-local effects arise due to the hadronic substructure of the photon. Because 
photon conversion into light partons is a genuinely long-distance process, the decay B — » X s 7 
is not a truly inclusive process, for which an expansion in local operators would apply. Indeed, 
from a conceptual point of view, it is as complicated as the semi-inclusive decay B — > X s h, 
with h denoting a specific light hadron. No analogous effects arise in semileptonic processes, 
since the conversion of heavy W bosons into light partons is a short-distance process. 

Effective field theories, such as soft-collinear and heavy-quark effective theory, provide the 
necessary tools to analyze inclusive B decays into light partons in the kinematical region of 
low hadronic invariant mass and large recoil energy, in which the hadronic final state is made 
up of a jet of collinear partons. For B — > X s ^ this is the endpoint region, where the photon 
has large energy E 1 m.b/2 in the B- meson rest frame. Effective field theories are systematic, 
taking into account all possible contributions to a given decay amplitude and describing them 
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in terms of well-defined, field-theoretic objects. This is especially important for radiative B 
decays, where the diagrammatic approach used in the previous decade has missed the largest 
source of non-perturbative uncertainty [39]. 

In this paper, we have shown that the B — > Xg'j photon spectrum in the endpoint region 
obeys the novel factorization formula (3). The first term in this formula has the structure 
familiar from semileptonic B decays. At each order in the l/m b expansion, it features products 
of hard functions Hi and jet functions Jj convoluted with soft functions Si. We refer to this 
term as the direct photon contribution, since the photon couples directly to the weak vertex 
in a local interaction. The two remaining terms in the factorization formula describe resolved 
photon contributions, in which the photon couples indirectly to the weak vertex via conversion 
into light partons. The partonic substructure of the photon is described in terms of a new class 
of jet functions Jj. These new terms appear first at order l/mb and arise from the contribution 
of operators other than in the effective Hamiltonian. The new soft functions Si entering 
the resolved photon terms contain non-localities in two light-cone directions. Only one non- 
locality is removed when the photon spectrum is integrated over energy to obtain the total 
decay rate. 8 As a result, we find that even the total decay rate receives non-local corrections of 
order AqcdA 7 ^- The new jet functions Jj, which are defined in terms of propagators dressed 
by Wilson lines, are complex quantities carrying calculable, perturbative strong-interaction 
phases. The soft functions, on the other hand, were shown to be real in by the use of heavy- 
quark symmetry and the invariance of the strong interaction under parity and time reversal. 
The impact of the new strong phases on CP violation in B — > X s ^y decay will be considered 
in a future publication. 

Phenomenologically the most important operators in the effective weak Hamiltonian are 
Qsg, and Q\. We have explicitly evaluated the l/mj, corrections to the B — > X s ^ photon 
spectrum arising from these operators, at tree-level in hard and hard-collinear interactions. 
This includes important contributions involving a hard-collinear gluon exchange, which carry 
a factor g 2 = 4na s . Our results are summarized in relations (12), which replace the relations 
(7) used in previous analyses of B — > X s 7 decay. The systematic methodology offered by 
the effective field-theory approach resolves a couple of puzzling features of the expressions 
(7), such as the appearance of the strange-quark mass in the expression for -F 8 p 8 art , or of large 
logarithms in the expressions for F p 7 art and -F 8 p 8 art . We point out that these features result from 
an improper separation of short- and long-distance physics. In our improved expressions (12), 
all long-distance physics is parameterized by well-defined hadronic matrix elements (the soft 
functions), while the logarithms entering the short-distance perturbative contributions contain 
(9(1) ratios of scales. At order l/m&, we find resolved photon contributions arising from the 
operator pairings Q$ g — Q$ g , Qi-y — Qs g , and Q\ — in the squared decay amplitude. We 
also prove that the resolved photon contributions arising from the Q\ — Q\ and Q\ — Q 8g 
operator pairings are suppressed by two powers of l/mb. Detailed analyses of the resolved 
photon contributions were presented in Section 4, which constitutes the main technical part 
of the paper. 

The non-perturbative soft functions, which are needed to describe the photon spectrum at 

8 As before, by "total" rate we mean the rate defined with a lower cut Eq on the photon energy that lies 
far outside the endpoint region, i.e., nib — 2E n >• Aqcd- 
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order l/m b in the heavy-quark expansion, introduce new sources of hadronic uncertainties in 
the description of the photon spectrum in the endpoint region. These functions can neither 
be extracted from experiment, nor can they be computed using lattice gauge theory (since 
they involve operators containing fields separated by light-like distances), and unfortunately 
they are not much restricted by constraints on their normalization and moments. Hence, 
there is a vast freedom in constructing phenomenological models for the soft functions, which 
often depend on several convolution variables. The resulting uncertainties will impact any 
extraction of \V u b\ via a combination of inclusive semileptonic and radiative decays. They will 
also affect the extraction of heavy-quark parameters such as m b , A, /i 2 etc. from moments of 
the B — > X s 7 photon spectrum. A dedicated analysis of the resulting uncertainties will be 
presented elsewhere. 

Our most important phenomenological result concerns the non-local power corrections to 
the B — > Xg'-f decay rate defined with a cut E 1 > E , where E is chosen to be far outside the 
endpoint region, mj — 2E ^> Aqcd- In this region the direct photon contributions reduce to 
local matrix elements, and deviations from the naive model of a free heavy-quark decay start 
at order l/m 2 c and 1 jm\ and are calculable in terms of well-known heavy-quark parameters. 
The resolved photon contributions, on the other hand, are still expressed in terms of non-local 
operators, whose matrix elements are of order Their contributions to the integrated 

rate can be parameterized in terms of three non-perturbative parameters, A17, A s ^ cc , and 
A 8 8, as shown in (103). In (104), these parameters are expressed in terms of convolutions of 
calculable jet functions with non-perturbative soft functions. Needless to say, it is very difficult 
to estimate the values of these hadronic parameters. Nevertheless, we have provided arguments 
suggesting that all three parameters are much smaller than the naive expectation ~ Aqcd- 
For the most important case of An, a detailed modeling of the corresponding soft function, 
taking into account the normalization conditions and moment relations we have derived in 
(49) and (53), suggests that A i7 is significantly smaller in magnitude than Aqcd, see (111). 
For the second- most important case of Ay| e c , we have provided two different arguments, based 
on the vacuum insertion approximation and on SU(3) flavor symmetry, suggesting that to 
a good approximation this hadronic parameter is proportional to the electric charge of the 
light spectator inside the B meson, see (113). While the parameter A^ x entering in this 
relation can indeed be of order Aqcd - see (114) and (116) - the weighting by the spectator 
charge reduces the corresponding contribution to the isospin-averaged decay rate by a factor 
(e M +erf)/2 = 1/6. Finally, as shown in (119), the parameter A 88 is suppressed by a charge factor 
e 2 s = 1/9, and its value is therefore bound to be much smaller than Aq C d- Our final estimates 
for the hadronic uncertainty from non-local 1 / rrib corrections in the theoretical prediction for 
the B — > X s 7 decay rate defined with the cut E 1 > 1.6 GeV has been given in (122) and (123). 
It depends on whether the contribution from Ayg CC is estimated using the vacuum insertion 
approximation or the current experimental value of the isospin asymmetry A -. We have 
emphasized that even a precision measurement of the isospin asymmetry would not help to 
reduce the uncertainty by a significant amount. Relation (124) shows that in this ideal case an 
irreducible uncertainty of about 4-5% remains. At present, we do not see any hope to reduce 
this error using well-controlled theoretical methods. 

The analysis presented in this paper applies without alteration, apart from some obvious 
substitutions, to the decay B — > X d ^. First, one needs to change the definition of \ q to V q bV* d 
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and replace m s in (34) and (35) by m d ps 0. Second, one has to replace the s-quark fields by 
ci-quark fields in the definitions of the soft functions /„ and /„ contributing to Fff F (E y , ji) in 
(11), and in the definitions of g^s and g^ contributing to F$(E y , fi) in (74). Notice, however, 
that none of these functions contribute to the integrated decay rate. 
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A Effective weak Hamiltonian 

We use the form of the effective weak Hamiltonian for B — > X s ^ decay as presented in [49], 



i.e. 



H e fi — —i= ^2 \ ( C\ Ql + C 2 Q\ + ^ Cj Qi + C 7j Q 7y + C 8g Q% g J , (A.l) 

v q=u,c ^ i=3,. ..,6 ' 

where \ q = V qb V* s , and the Wilson coefficients depend on the scale /i at which the operators 
are renormalized. The explicit form of the operator basis is 

Ql = (qb)v-A (sq)v-A , Ql = (qibj)v-A (sjq^v-A , 

Q3 = {sb) V -A ^2 1 i.QQ)v-A , Qi = {Sibj) V -A ^2 1 (QjQi)v-A , 

Q5 = {sb) V ~A ^2 1 (Vtiv+A , Q& = (Sibj)v-A ^2 1 iWi)v+A , 



(A.2) 



Qr, = -^r + 75)^ , Qs 9 = *v(l + 75)C^6 , 

where % and j are color indices, and for the penguin operators a summation over quark flavors 
q — u,d, s, c, b is implied. We use the short-hand notation (qb)v^A = <?7 M (1 =F 75)^ etc. Our 
sign convention is such that iD^ = id^ + gT a A a ^ + ee q A^, where T a are the SU(3) color 
generators, and e q are the quark electric charges in units of e. 
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B Input parameters 



Here we collect the input parameter values used in the numerical analysis in Section 7. As our 
default choice for the factorization scale \x entering the Wilson coefficients, the strong coupling 
constant, and the various hadronic quantities, we take the hard-collinear scale \i = 1.5 GeV, 
which is indeed a scale of order ^/ rrib Aqcd- This is an appropriate scale choice, given that we 
neglect RG evolution effects. 

The 6-quark mass enters our expressions either via the photon energy (E 1 ~ rrib/2 near the 
peak of the spectrum) or as the heavy-quark expansion parameter. It is therefore appropriate 
to adopt a low-scale subtracted heavy-quark mass, such as the mass defined in the shape- 
function scheme [29]. Specifically, we use nib = 4.65 GeV. The charm-quark mass enters as a 
running mass in charm-penguin diagrams with a soft gluon emission, which are characterized 
by a hard-collinear virtuality. We therefore use m c = m c (//) defined in the MS scheme, with 
\x = 1.5 GeV fixed as described above. This corresponds to the choice adopted in [5], and 
following these authors we set m c (//) = 1.131 GeV. Finally, for the strange-quark mass we 
take m s = rrib/50, which is the value commonly adopted in the literature on B — > X s 7 decay. 

We also need input values for some HQET matrix elements. The parameters A2 and p\ s 
are extracted from a global fit to B — > XJ v experimental data by the Heavy Flavor Averaging 
Group (HFAG) [7]. Unfortunately, in many cases these and other parameters are extracted 
from a combined fit to B — > X c lu and B — > X s ^, an approach that was criticized in [81]. 
Only recently HFAG has started quoting also values obtained using only semileptonic data. 
The most recent results are A 2 = (0.12 ± 0.02) GeV 2 and p 3 LS = (-0.17 ± 0.09) GeV 3 [82]. For 
simplicity, we always use the central values for these quantities. For the first inverse moment of 
the _B-meson light-cone distribution amplitude, we take the range 250 MeV < \ B < 750 MeV, 
which covers predictions obtained using QCD sum rules and other methods [72-77]. Finally, 
to the level of accuracy of our calculations, the parameter F can be extracted from the relation 
F = f B ^/MB, and using f B = (193 ± 10) MeV [83] we obtain 0.177 GeV 3 < F 2 < 0.217 GeV 3 . 

C NNLO matching of to SCET 

In this appendix we present the matching of the effective weak Hamiltonian operators Q77, 
Qg g , and Q c { u onto SCET up to NNLO in the expansion parameter V% with A ~ Aq CD /m fe . 
Although there is a large number of possible operators, only some of them are needed in 
practice. One subset, which was presented already in Section 4.1, is needed for the study 
of the resolved photon contributions at tree level. Another subset is needed for the analysis 
of the power corrections to the direct photon contributions at 0(a s ). In the first part of 
this Appendix we perform the matching at tree level. In the second part we include also the 
contribution of one-loop quantum fluctuations. 

C.l Tree- level matching 

We begin with the current-type operators and Q$ g . At LO, with a scaling of A 5 / 2 , the 
operator is the only one which gives a contribution to B — > X s ^. The contribution of 
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the operator Q 8g begins at NLO, i.e. C(A 3 ). We then perform the tree-level matching of Q c { u , 
whose contribution begins at the NNLO, 0(\ 7 / 2 ). For simplicity, we denote LO, NLO, and 
NNLO operators by superscripts (0), (1), and (2), respectively. 



C.l.l Matching of Q 7l 

The operator Q 7l in the weak Hamiltonian is given by 

Qr 7 = -fP^(l + = -^s (1 + 7s)& , (C.I) 

where it is assumed that the photon is real, i.e., it is transversely polarized. The tree- level 
matching of Q77 can be read off from [43]. As was done there, we separate Q77 into "A" and 
"B" terms, according to whether they contain hard-collinear gluon fields or not. Suppressing 
the — f^r e~ tnibV ' x factor, Q 7l is matched onto the operators 



(0) 

77A 


, f 
S/ic 2 


(1) 

77A 


S/ic 2 


(2) 
77A 


S/ic 2 



n • x x^ x\ ^ „ , ill) , 

+ &o (1 + 75) xf x D„h + f te I [ih4T\ (1 + 75) A , 

1 in- d 



(C.2) 



and 



QZ = he [ - ' d4T] 9 Ac±(i + 75)^, 

171 b 

Q?]b = -Lc I [in ■ d4T] (1 + 75) 9 n ■ A hc h + £ hc [ — ' g n ■ A hc {\ - 75)^ 



+ Lc [ — ' HT] 94hc±(l + l,)x^D,h (C.3) 

171 b 

-^-[m.fl^Kl + Ti)^ ft 

- 4c ^4^c± — ^= [^±41 m ] (1 + 7s)/i , 

^ in - d 

where the covariant derivative -D M only contains soft fields. Note that the hard-collinear fields 
are not sterile (see the discussion in Section 3), so they still couple to soft gluons. 



53 



C.1.2 Matching of Q 8g 



The operator Q 8g in the weak Hamiltonian is given by 



' 8tt 2 



^(1 + 75)^6. 



(C.4) 



We can discard the non-abelian part of G^", since we only work to first order in g. A priori, 
we can match the s quark onto either a hard-collinear or anti-hard-collinear quark and the 
gluon onto a hard-collinear or anti-hard-collinear gluon. We do not consider matching the s 
quark onto a soft particle, since it will only give rise to a contribution beyond NNLO. Also, 
the gluon and the s quark cannot both be anti-hard-collinear, since the necessary conversions 
would lead to a suppression beyond NNLO. The three remaining cases will be considered in 
turn. 



SCET 



he 

| 

» he 



•-SCET 



+ 



he 



s s 
he" 



WW 



v/VWV 



f 



■9- 



Figure 10: Graphical illustration of the tree-level matching procedure for a Q 8g contribution, 
showing how a resolved photon contribution arises after a SCET Lagrangian insertion. 



For the case of an anti-hard-collinear gluon and hard-collinear s quark in Q Sg , the conver- 
sion of the anti-hard-collinear gluon onto a photon is 0(A) in power counting. This process is 
illustrated in Figure 10. In this case Q$ g is matched onto a single operator (again suppressing 
the overall -a^k e -im b vx f actor ). 

QsgM gluon = ^ c f ^ n ' WhZil ( X + , followed by 0(A) conversion. (C.5) 

This is the first example of an operator that gives rise to a resolved photon contribution. The 
conversion is displaced along the light cone by the anti-hard-collinear propagator. 

For the case of a hard-collinear gluon in Q 8g , we have to distinguish between the different 
scaling of the components of . We will need the following components of G^ v : 

0(\ 1/2 ) : 2 | [in- 04 J 

0(A) : 2 [i$ ± 4 ± - id ± ■ A ± ] + (|| - |0 [in • On • A] 

If the s quark is matched onto a hard-collinear field, we need to have the anti-hard-collinear 
photon emitted from the b or s quark lines, as shown in Figure 11. This would match onto 
SCET operators that contain both a hard-collinear gluon and an anti-hard-collinear photon 
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SCET 




Figure 11: Graphical illustration of the tree- level matching procedure for the operator Q Sg for 
the case when the s quark is matched onto a hard-collinear field. 



suppressed by mj, for the 6-quark line, and by n ■ d for the s-quark line. We also need to 
consider the three possibilities for and the multipole expansion of the heavy-quark field. 
In total we have 

QShcquark = & C — ^=^47 [ifi • dJfi hcA ] (1 + lb )h , 

in ■ a 



(2) 

t> 8g, he quark 



c 



i 



e d e4l m [m • dJfi hcA ] (1 + ^x^D^h 



hc~ ^= 

in • d 

- £ hc ^^J- [in • a^eJe^Hl + 75)^ 
in • d ^ m b 

1 $ 

+ Che— ^e d e41 m - [i<^^ c± - id ± • .A^J (1 + 75) /i 
in - d 1 



(C.7) 



- Che [i$±4he± - i<9± • A hcA ] - — e d e47(l + ^ 5 )h 

2m b 



+ Che 
+ Che 



*=e d e4T^ 

in - a z 



-n- dn- A hc 



(1 + 75)/* 



2m h 



-n ■ dn ■ Ahc 



e d e41 m (l + 7 5)/i. 



If the s quark is matched onto an anti-hard-collinear field, it can only be converted to an 
anti-hard-collinear photon and a soft s quark. The conversion costs us A 1 / 2 , so the lowest-order 
operator possible is 0(\ 3 ). Considering all the possible structures for G^ u and the multipole 
expansion, we find 



Q (1) 
Q 



8g, he quark 
(2)_ 

8g, he quark 



V n ' d4-hc±\ (1 + 75)^ , followed by 0(X 1 ^ 2 ) conversion, 

ilTe^ \™ ■ d4he±] (1 + ^x^Dph , followed by £>(A 1/2 ) conversion, (C.8) 
+ £hc~ [i$±4hc± — id± ■ A hc ±] (1 + 7 5 )/i , followed by 0(X 1 ^ 2 ) conversion, 

l + 7 5 )/i, followed by £>(A 1/2 ) conversion. 



ihc 



i 
2 



-n-dn- A hc 
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C.1.3 Matching of Q\ 

To simplify the notation we write the operator as Q\ = sTiqqT 2 b, where Ti (g) T 2 = 7 M (1 — 
75) <8> 7^(1 — 75)- At tree level, the light quarks can only be matched onto hard-collinear or 
anti-hard-collinear fields. A matching of any of the light quarks onto even one soft field would 
lead to a suppression of 0(X 4 ). As a result, Q\ is matched at NLO, before taking into account 
any conversions. When there is more than one anti-hard-collinear quark field, no conversion 
is allowed at tree-level. Hence, we are left with the following cases (suppressing the e~ tmb v ' x 
factor). 



s = £hc> q = £ho q = £hc- We need to consider an attachment of one anti-hard-collinear 
photon emitted either from the heavy or one of the hard-collinear quark lines. This leads to 
four possible 0(A 7 / 2 ) operators: 

Zm b 1 in- d (C 9) 

+ IhcTMte e d e4T7;—^=r 2 ^ hc - ^ hc T 1 h^ hc T 2 ^ e d e41 m — L^ c ■ 
1 in- d 1 in - d 



s = €hc q = €hE, q = ihc and s = £ hc , q = £ hc , q = In this case we need to 

convert the anti-hard-collinear quark to a photon via £^ — > A" 11 + q and ^ — > A c ™ + q, which 
is 0(X 1 ^ 2 ). Therefore, we can only have 

^ c r 1 /i^r 2 Ck , followed by C(A 1/2 ) conversion, 

(CIO) 

ZhcTihthcTzZhZ, followed by 0(X 1 ) conversion. 



s = £^r, q = ^hc, q = £ hc : Supplemented with SCET Lagrangian for ^ — > A™ + q, only 
one NNLO operator is possible: 

lhFihl hc T 2 £ hc , followed by C(A 1/2 ) conversion. (C.ll) 



C.2 Loop matching 

We now perform the matching including the contributions of loops. These are only relevant for 
Qi, where the two up-type quarks are contracted and a number of gauge bosons are emitted 
from the internal lines. The contribution of three gauge bosons would lead to further power 
or loop suppression, so we only need to consider one or two bosons. 

The one gauge boson loops are easy to analyze. In the NDR scheme only Q 5 and Qq give 
a non zero contribution. Furthermore, this contribution only modifies the coefficients Q77 and 
Q 8g to Cf [84], with 

6 6 

c?7 — + ^2 > = c 8g + ^2zid , (c.12) 

i=i i=\ 



56 



where y^ Zi depend on the scheme. They vanish in the 't Hooft-Veltman scheme, while in the 
NDR scheme the non zero ones are 2/5 = —1/3, y§ = — 1, z^ = 1. The effective coefficients are 
regularization-scheme invariant. The contribution of the one-boson loop would therefore be 
to change CV 7j8s to C^ 8g . 

A more involved contribution arises for the loops with two external bosons, which is the 
main focus in this subsection. These can only be one photon and one gluon. Two external 
gluons would lead to further power or loop suppression. 

As usual, the b quark is matched onto the heavy quark field h. Since there is already 
a photon in the operator, the s quark cannot be anti-hard-collinear. On the other hand, a 
soft s quark would lead to power suppression, so it must be hard-collinear. There are two 
possibilities for the gluon emitted from the loop: it can either be hard-collinear or it can be 
soft. If the gluon is hard-collinear, the loop momentum is hard. If the gluon is soft, the loop 
momentum is anti-hard-collinear. For the first case, a photon and a hard-collinear gluon, we 
need to calculate the loop diagram in QCD in order to perform the matching. For the second 
case, a photon and and a soft gluon, one would need the tree-level matching of Qi onto the 
operator of (25). The conversion of the two anti-hard collinear quarks to a photon and a 
soft gluon would be calculated in SCET. Alternatively, one can calculate the process in QCD 
and use the fact that the two calculations are equivalent, since only one momentum region, 
anti-hard-collinear, contributes in this case. 

We have explicitly calculated the appropriate QCD one-loop diagram, arising from the 
four-quark operator (siY 2 qj) ((^Ti^)- Alternatively, the result can be read off from the two 
gluon calculation in [85] , where one of the gluons is replaced by a photon. Using the notation 
of [85] the amplitude is given by 



= ee q g 

where 



A = ^^(T a ) mn s m r 2 A^r 1 b n er(q 1 )e* au (q2)S ij 5 kl , (C.13) 

47T 47T 



V = [(4r - l)F(r) - 4r] 



2r m„ 



+ [F(r) - 1] i W {q[ - q p 2 ) 7V 



2r F(r) 
H [F(r) - 1] [q^ie paav q p q 2 - qviep^qlql] 7°7 5 «W^i<?27 5 

171^ TTlg 



(C.14) 



Here m q and e q are the mass and charge of the quark in the loop, qi,€\ {q 2 ,^2) are the 
momentum and polarization of the photon (gluon), r = m 2 q /q 2 — ie, q 2 = (qi + q 2 ) 2 , and F(r) 
is the penguin function defined in (9). Alternatively, one could write in a gauge- invariant 
notation 



.VrI"(</iK"'(</2) = [CI - ! )/•'(/■) [ ^^ F ^ G U + ^ G ^ FP ^ 

9 q (C15) 

+ [1 - F(r)} 1 (G^Fi* + F, a G a ^) iq a lpl , . 

Here we are using the convention 

2 
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and the fact that for an external gluon q 2 ■ e*£ = 0. 



C.2.1 Matching with A ^ and Af uon 

The diagram with an anti-hard-collinear photon and a soft gluon emitted from the internal 
line already contributes at 0(\ 7 ^ 2 ). Since A s gluon scales homogeneously ~ (A, A, A), it is natural 
to use the gauge invariant form for the resulting operators, rather than decomposing the gauge 
fields into their light-cone components. Furthermore, only the axial part of A^ v in (C.14) yields 
a non-zero result, as the Tj in (C.13) are the usual V — A Dirac structures, when matching 
Qi. Therefore we only need 



Ve^Oer (? 2 ) = - (g%F^ + F, a G a ^) iq a ^ 5 [l - F(r)} 

2 (C - 17) 



which follows from the fact, that q a F fia vanishes at the lowest order in A. 

For Qi the loop can consist of any up-type quark that is not integrated out in the weak 
effective Lagrangian. When matching onto SCET the u quark should be taken to be massless, 
so we can replace F(m 2 /q 2 ) by 0. For charmed quarks we have m 2 ~ 771{,Aq CD , which is of the 
same order as q 2 . As a result, F(m 2 /q 2 ) should not be expanded for c quarks. 

In position space, we find that at NNLO the Q\ operators are matched onto 

Q? (2) (x) = (0) TV(1 - l,)Kx)e- im ^ 

Qf\ x ) = (0) Ux) TV(1 - l B )h(y)e-**™ J ^ (^ ^^ ( C ' 18 ) 



1 

x 



(<7i + Q2) 2 + it 



2 

- re 



{fi + CL^gGlMe^FM, 



where we show explicitly the dependence of the momentum-space F pa and G a m on q\ and q 2 - 



C.2.2 Matching with and Afl™ 

In this case q 2 is hard, i.e. q 2 ~ m 2 . Therefore F(m 2 /q 2 ) can be expanded around zero for 
charm quarks as well as for up quarks. The first order correction resulting from this expansion 
gives a power suppressed contribution, so we can just set F(m 2 /q 2 ) to zero. Depending on the 
polarization of the gluon field we can get either an NLO or an NNLO operator. We need to 
also include corrections from the multipole expansion of the heavy-quark field and subleading 
matching on the hard-collinear quark. In total we find that Q\ is matched onto (suppressing 
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an overall factor £Sg e -i"*6 

Qt (1) = e ± , u [Alc± n ■ dA? 1 (1 - 75)^ , 



2 

g? (2) = skccw^r • ^fccj.] 7^(1 - 75)/» 

+ \Zhce± l u,[ A T v n ■ dn ■ A hc ] 7^(1 - 75)^ (aig) 
in ■ 

+ a c f Ci^I^n • BAT 1 (1 - TfeX^ » 



where ej.^ = \e a p lll/ n ai nP . 
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